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Abstract. Abstract. In this work we apply the Poincare-Cartan formalism of 
the Classical Field Theory to study the systems of balance equations (balance 
systems). We introduce the partial k-jet bundles Jp(ir) of the configurational 
bundle it : Y — > X and study their basic properties: partial Cartan structure, 
prolongation of vector fields, etc. A constitutive relation C of a balance system 
BC is realized as a mapping between a (partial) k-jet bundle Jp(^r) and the 

extended dual bundle A^t^ n ~^^Y similar to the Legendre mapping of the 
Lagrangian Field Theory. 

Invariant (variational) form of the balance system Be corresponding to a 
constitutive relation C is studied. Special cases of balance systems -Lagrangian 
systems of order 1 with arbitrary sources and RET (Rational Extended Ther- 
modynamics) systems are characterized in geometrical terms. Action of auto- 
morphisms of the bundle tt on the constitutive mappings C is studied and it 
is shown that the symmetry group Sym(C) of C acts on the sheaf of solutions 
Sole of balance system Be- Suitable version of Noether Theorem for an action 
of a symmetry group is presented together with the special forms for semi- 
Lagrangian and RET balance systems and examples of energy momentum and 
gauge symmetries balance laws. 
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1. Introduction. 

Systems of balance equations (balance system, shortly BS) for the fields j/ M 
accompanied by the proper constitutive relations C are the basic tools of Continuum 
Thermodynamics. In this work we develop a variational (in the sense of Poincare- 
Cartan form) theory of balance systems. 

In Sec. 2 we introduce basic notions: k-jet bundles, contact decomposition of 
forms, Poincare-Cartan formalism of order 1, balance systems. In Sec. 3 the partial 
jet bundles Jp (tt) of configurational bundles tt : Y n+m — ► X n are introduced as the 
appropriate domains of the constitutive relations of the corresponding field theory. 
These bundles corresponds to an almost product structure on the base manifold 
M (space-time decomposition is the main example) . Partial Cartan structures and 
prolongation of vector fields from Y to the partial jet bundles are studied. In Sec. 4 
the constitutive relations C are defined as the mappings C : Jp (tt) — > Ay^ n+1 ^Y 

to the bundle AoJ^~ F of n + (n + l)-forms on Y annulated by two 7r-vertical 
tangent vectors factorized by the bundles of similar forms annulated by one vertical 
tangent vector. Lifting C : J^tt) — * A2 + ^ n+1 ^Y of such a mapping induces the 
n + (n + l)-form on the partial jet bundle - Poincare-Cartan form of consid- 
ered field theory. Special forms and examples of such CR are presented: lifting 
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covering constitutive relations, semi-Lagrangian, RET, relations defined by the La- 
grangian and a Dissipative potential, etc. In Sec. 5 the invariant form of a balance 
system corresponding to a constitutive relation C is introduced and studied. Cases 
of Lagrangian field theories of the first order and of the Rational Extended Ther- 
modynamics as a specific case of the balance systems are described. We study the 
admissible variations of the fields yfi splitting the invariant equation into separate 
balance equations. An appropriate modification of the source term allows us to 
remove conditions on the admissible variations. 

Action of automorphisms of configurational bundle n : Y — > X on the constitu- 
tive relation C is defined and its effect on the corresponding Poincare-Cartan form 
and on the sheaf of solutions of the balance system is studies is Sec. 6. 

In Sec. 7 we prove the appropriate version of the first Noether Theorem associat- 
ing the new balance laws with an infinitesimal symmetries of a constitutive relation 
C and determine when such a balance law is actually, the conservation law. Exam- 
ples of the energy-momentum balance law and that of the gauge symmetries are 
presented. In Sec. 8 we sketch the application of present scheme to the search and 
classification of the "secondary balance laws of a given BS, including the entropy 
balance law. As an example, we present classification of secondary balance laws for 
the Cattaneo heat propagation balance system and the constitutive restrictions on 
the CR C that follows from the II law of thermodynamics. 

Short exposition of main results of this paper was presented in the Conference 
Proceedings [27] . 



2. Settings and the framework 

2.1. Notations and preliminaries. Throughout this paper it : Y — > X will be a 
(configurational) fibred bundle with a n-dim connected paracompact smooth (C 00 ) 
manifold X as the base and a total space Y, dim(Y) = n + m. Fiber of the bundle 
7r is a m-dim connected smooth manifold U. 

Base manifold X is endowed with a (pseudo)-Ricmannian metric G. Volume 
form of metric G will be denoted by n. In this paper we will not be dealing with 
the boundary of a base manifold X, in applications X can be considered as an open 
subset of R n or as a compact manifold. As an basic example, we consider the case 
where X = T x B is the product of time axis T and an open material manifold (or 
a domain in the physical space) B. 

We will be using fibred charts (W, x % , y M ) in the bundle tt. Here (it(W),x 1 ) is a 
chart in X and y 11 are coordinates along the fibers. Tangent frame corresponding to 
the local chart (W, x l , y^) will be denoted by {pi — d x % , = d y ^) (shorter notation 
will be used in more cumbersome calculations), corresponding coframe - (dx l , dy^). 

Introduce the contracted forms r/i — ig i n, r]ij = ig j ig i r\. Below we will be 
using following relations for the forms T]j, rjji (here and below Xq = ln(\/\G\)): 

{dx j A rn — 5jn, 
dr/i = Ag,<?7, 
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Sections s : V — > Y, V C X oi the bundle tt represent the collection of (classical) 
fields defined in the domain V C X. Usually these fields are components of 
some tensor fields or tensor densities fields (sections of "natural bundles", [1]). 

For a manifold M we will denote by tm : T(M) — » M the tangent bundle of a 
manifold M, by V(M) C T(M) the subbundle of of vertical vectors, i.e. vectors 
£ G T y (M) such that T M * v i = 0. 

Denote by A r (M) the bundle of exterior r-forms on the manifold M and by 
(A* = ©A r (M), d) the differential algebra of exterior forms on the manifold M. 

2.2. The k-jet bundles J fe (7r). Given a fiber bundle tt : Y — > X denote by J fe (7r) 
the k-jet bundle of sections of the bundle w, [12]. Denote by iTkr ■ J {n) - * 
J r (7r), k ^ r ^ the natural projections between the jet bundles of different order 
and by ir k : J k — > X the projection to the base manifold X. Projection mappings 
■^feffe-i) : J k { 7r ) ~* J (ft) in the tower of k-jet bundles 

. . . -» j fe w -> J fc_1 (7r) Y" —► 

are affine bundles modeled by the vector bundle /\ k T*(X) ®^jk-i^ V(ir) — * 
J k - l {ir). 

Denote by </°°(7r) the infinite jet bundle of bundle tt - inverse limit of the pro- 
jective sequence KMk-x) '■ J k ( n ) ~ * J k ~ 1 ( 7r )- Space J°°{-k) is endowed with the 
structure of inverse limit of differentiable manifolds with the natural sheaves of 
vector fields, differential forms etc. making the projections ir^k : J°°(7r) — > J k (ir) 
smooth surjections. See [12l [28] for more about structure and properties of k-jet 
bundles. 

For a mutliindex / = {ii, . . . , i n }, £ N denote by d 1 the differential operator in 
C°°(X) d 1 f = d l x \ ■ . . . 9*" . To every fibred chart (V, x l 7 y^) in Y there corresponds 
the fibred chart {x\y»,z^ |/| = £ a *« ^ k) in the domain V k = n^(V) C J k (n). 
This chart is defined by the condition Zj(j k s) = d 1 s ,i (x). 

For k — 1, . . . , oc the space of k-jet bundle J k {^) — > X is endowed with the 
Cartan distribution Ca k defined by the basic contact forms 

uj» = dy^ - z?dx\ ...w% = dzf- zf +lj dx j , |/| < k; dz?, \I\ = k. (2.2) 

These forms generate the contact ideal C k C A*(J fe (-7r)) in the algebra of all exte- 
rior forms. Denote by I(C k ) the differential ideal of contact forms. This ideal is 
generated by the basic contact forms (2.2) and by the forms dzj , |/| = k, see [11] . 
In the case where k = oo the basic contact forms dzj are absent form the list of 
generators of ideal I(C k ). 

A p-form is called /-contact if it belongs to the 1-th degree of this ideal (C k ) 1 C 
A*(J fc (7r)). 0-contact forms are also called horizontal (or 7Tfc-horizontal ) forms 
(or, sometimes, semi-basic forms). We denote by kCon the fc-contact forms that 
appear in calculations. For k = 1 we will omit index 1. 

Let ^ s < k. A form v G A*(J fc (-7r)) is called 7Tk s -horizontal if it belongs to 
the subalgebra C*°°( J fc (7r))7r* s A*( J s (vr)) C A*(J fc (?r)). 

Remind now the following basic result (D. Krupka,[9]) 

Theorem 1. Let aformv G A 9 (J fe (7r)) be ir^k-i)' horizontal (i.e. v = Tt'klk-i) 1 '*, v* G 
A 9 (J fe - 1 (7r)). Then there is unique contact decomposition of the form v 



V = Vq + V\ + . . . + V, 



(2.3) 
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where z/j,0 ^ i ^ q is a i-contact form on J k (n). Form vq is called the horizontal 
part of the form (and of the form v as well). 

Related to the contact decomposition is the decomposition of differential operator 
d as the sum of horizontal and vertical differentials dh,d v : for a q-form v G 
A*(J fe )(7r)) its differential dv lifted into the J fc+1 (7r) is presented as the sum of 
horizontal and contact (vertical) terms: 

dv = d^v + d v v. 

,see [5J[T3]. Operators dh,d v are more naturally defined in the space A*(J°°)(7r)). 
We recall that these operators have the following homology properties 

dl = dy = d v d h + d h d v = 0. 

In particular, for a function / G C°° (J°° (%)) (depending on the jet variables Zj 
up to some degree, say \I\ ^ k), 

df = (dj)dx* + ( 2 -4) 

where 

dif = d x if+ tf+ufyf ( 2 -5) 

is the total derivative of the function / by x\ The series in the formulas (2.4-5) 
contains finite number of terms: |7| ^ k. 

2.3. Lagrangian Poincare-Cartan formalism, k=l. Here we remind the basic 
notions of infinitesimal form of Lagrangian Field Theory of the first order (on J 1 (ir)) 
based on the use of Poincare-Cartan form,[H [51 [T5] . 

Volume form rj permits to construct the vertical endomorphism 

S v = (dy»-z?dx i )Ar ]j ®-^ (2.6) 

which is a tensor field of type (l,n) on the 1-jet bundle J 1 (7r). 

For a Lagrangian n-form Lrj, L being a (smooth) function on the manifold J 1 {it) 
the Poincare-Cartan n-form are defined as follows: 

e L = Lr ) + S;(dL), (2.7) 

where S* is the adjoint operator of S v . In coordinates we have 

®L = (L-Zi-^)v+-Qjxdy»A m , (2.8) 

An extremal of L is a section s of the bundle 7r such that for any vector field £ on 
the manifold J 1 (7r), 

(j 1 ( s )y(i i de L ) = o, (2.9) 

where j (s) is the first jet prolongation of s. A section s is an extremal of L if 
and only if it satisfies to the system of Euler-Lagrange Equations (see, for instance, 

E2U5J) ' 
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2.4. Bundles K r p Y and canonical forms. Introduce the bundle AL(Y) of the 
exterior forms on Y which are annulated if p of its arguments are vertical ( [101116] : 

/eA;(r)« l6 .., £ / = 0, ti€V(Y). 

We will be using these bundles for r = n, n + 1, n + 2 and p — 1,2. 

Elements of the space A™F are semibasic n-forms locally expressed as p(x, y)rj. 
Elements of the space A%Y have, in a fiber coordinates (x l , y^) the form 

p(x,y)r)+p] Ji dy IJ Ar] t . 

This introduces coordinates (a; 1 ,?/^,^) on the manifold A™y and (x z , 
on the manifold A^y 

For the case where r = n + 1 the forms dy^ A 77 form the basis of fibers of the 
bundle A' 2 l+1 (r) -> F while the bundle A" +1 (y) is zero bundle. 

Introduce the notation 

A»+(™+ l )(y) = A™(y)©A™+ 1 (y),p = i,2, 

for the direct sum of the bundles on the right side. 

It is clear that Af (y) C A|(y). Therefore we have the embedding of subbundlcs 
A" +(n+1) (y) C A" +(n+1) (y) and can define the factor-bundle 
A „+(„+i) (y) = A «+(«+i)( r ) /A «+(«+i)(y) j 

with the projection q : A£ +(n+1) (Y) -> A£+ (n+1) (y). 

On the bundles A^ + + 00 there is defined the canonical form ([H]) with the 
coordinate expression 

e«+(«+ 1 ) = e « + e «+i = a % + P/t d tf " a r,. (2.11) 

On the factor-bundle A^J^ n+1 \Y) the form 

e; l + (n+1) = p'.dy" A m + p fi dy" A Ti 

is defined mod r\. 

2.5. Legendre transformations, k = 1. Let i G C°°(J 1 (7r)) be a Lagrangian 
function. We define the fiber mapping (over Y) legL '■ J 1 ^) — * A^y, as follows: 

Zeffi(jia))(Jr x ,...,X B ) = (e^-x^X!,...,!^, 

where j^s G J 1 (tt), X G T S ^ X \Y and X, G Tji s t x \J (w) are such that 7r*(JQ) = Xj. 
In local coordinates, we have 

Ze 5L (x\?/\zf) = (x i ,y'*,p = L-z!t-^ jX ,pl = -^). 

The Legendre transformation Legr, : J 1 (7r) — > A™^ l 1 < ' n+1 ' ) (y) is defined as the com- 
position LegL = Q legL- In coordinates, 

BT 

Leg L (x\y»,z?) = (x\y»,pi = — T ). 

Recall [16] that the Legendre transformation Legz, : J 1 (tt) — » A^f (y) is a local 
diifeomorphism if and only if L is regular, i.e. when the vertical Hessian d ^g z „ is 
nondegenerate. 
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2.6. Con-differential. Below we will be using the con-differential d used in topol- 
ogy, see. for instance, [7]. 

Let M n be an n-dimensional smooth manifold. Con-differential is defined on 
couples of exterior forms a k + f3 k+1 e k k+l - k+1 \M): 

j. A*+(*+i)(M) = A k (M) © A k+1 (M) -» A(' £ + 1 )+( fe + 2 )(M) = A k+1 {M) © A fe+2 (M) : 
J(a fe + /3 fe+1 ) = ((-da + p) + d/3). 

(2.12) 
Lemma 1. rforf = 0. 
Proof. We have 

JJ(a fc + /3 fc+1 ) = J((-da+/?) + d/3) = [-d(-da+#) + d0\ +d{d/3) = -d(3 + d(3 + 0. 

□ 

Some other properties of this differential and of the corresponding complex 
( A fc+(fc+i)(M) = A k {M) © A k+1 (M),d) are presented in [25]. 

2.7. Balance systems. Let the base manifold X be the material or physical space- 
time X with (local) coordinates x 1 — t, x ,A — 1, 2, 3. Typical system of balance 
equations for the fields is determined by a choice of flux fields i* 1 ', /i = 0, 1, 2, 3 
(including the densities for i = 0) and source fields 11^ as functions on a jet space 
J fe (7r). Choice of the functions F^II^ s C°°(J fc (7r)) is codified in physics as the 
choice of the constitutive relation of a given material, [22] [23]. After such a 
choice has been done, the closed system of balance equations for the fields y % 

(F; o j k s).^ = n M o j fc s , M = 1, . . . , m (2.13) 

can be solved for a section s : V — > Y if one add to the balance system the 
appropriate boundary (including initial) conditions. 

Introducing the horizontal forms in J k {n) - = FJ^rji, = ILj77 these equations 
can be written in the form 

i 1 *( s )[dF M -ry = 0. 

Using the con-differential d (see above) the balance system (2.13) can be written 
in the compact form 

j 1 *(s)d[F tt + U tt }=0, i = l,...,m. (2.14) 

Remark 1. Similarly to the definition of a conservation law on a contact manifold 
given in [3] one can define a balance law of order k as a n + (n + l)-form a = 
F n + n™ +1 on the space J fc (7r) (which usually is 7Tfc-horizontal) where F n is defined 
mod closed forms from Z"(J fc (7r)) such that 

da = dF n - U n+1 G I{C k ). (2.15) 

In Continuum Thermodynamics balance laws typically are present as the closed 
system of equations for the dynamical fields. It is shown below that the Poincare- 
Cartan formalism allows to generate the system of balance laws in the number equal 
to the number of dynamical fields y^. 
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Example 1. Five fields fluid thermodynamical system. In the 5 fields fluid 
thermodynamical system ([221I23]), the velocity vector field v, density scalar field p 
and scalar field of internal energy e (or total energy e=^p|M| 2 + £) are considered 
as the basic dynamical fields y M . Corresponding balance system has the form of the 
five balance equations - mass conservation law, linear momentum balance law and 
the energy balance law: 

d tP + d x n(pv B ) = 0, 

d t {pv A ) + d x B(t p v A v B + t AB ) = } A , A = 1,2,3, (2.16) 
dt{hp\\v\\ 2 + e) + x b ((±HM| 2 + e)v B + t B v c - h B ) = f B v B + r. 

Here t is the (l,l)-stress tensor, f is the 1-form of the bulk force, h is the heat flux 
vector filed and the scalar function r - radiation heating source. 

Example 2. 



3. Partial k-jet bundles 

In this section we introduce the "partial jet bundles" of the configurational bun- 
dle 7r : Y — > X - factor bundles of the conventional jet bundles J k {^). These 
bundles are convenient for the description of balance systems whose constitutive 
relations depends on some but not all derivatives of fields y^. We start with the 
definition of the partial 1-jet bundle defined by a subbundle K C T(X). 

3.1. Bundle Jj^tt). 

Definition 1. Let ttxy ■ Y — > X be a fiber bundle and let K C T(X) be a 
subbundle. 

(1) Let x G U C X, s 1 ,s 2 G T(U,tt)\si(x) — s 2 (x). Sections s-y,s<2 are called 
K -equivalent of order 1 at a point x G X : Si ~k x s 2 if Si* x \k x = 

S2*x\K x - 

(2) J\r x (tt)- space of classes of ~k x at a point x G X. 

(3) Jhi 71 ) = ^xexJx x {ir) - the space of partial 1-jets of sections ofn. 

In the next proposition we collect basic properties of bundles Jh (tt) . Proof of 
this Proposition is straightforward ( see [2"5]). 

Proposition 1. (1) Bundle Jk-(tt) Y is the affine bundle modeled on the 
vector bundle n*(K*) ® V(w) — > Y. 

(2) There is a canonical surjection of affine bundles wk '■ J 1 (it) — > J^-(7r). 

(3) LetT(X) = K(X)®K'(X) be an almost product structure (AP), then there 
is the commutative diagram 

ji W jI kM 

W K'\ W10 k\ 

t1 I \ "10 K ' v 

which realizes the conventional 1-jet bundle J 1 (it) as the fiber product of 
the bundles Jk(tt), Jici 71 ')- 
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(4) Let T(X) = K © K' be an integrable AP-structure, (x^,x k ) - local chart 
integrating the AP — structure ( i. e. K x is the linear span of vector fields 
dj, K' - is the linear span of vector fields dk). One forms 

oj" = dy» - J] zfdx j 
j 

generate the "partial" Cartan distribution Ck on J^-(7r) in the sense that a 
section q of the bundle J^(tt) is the (partial) 1-jet of a section s : X — > Y 
if and only if q* (w^Ik = 0, V/i = l,...,m. 

(5) All constructions above are covariant with respect to the automorphisms 4> 
of the bundle n such that projection <f> of automorphism cj> to the base X 
leaves the AP structure T — K K' invariant. 

Remark 2. Partial Cartan structure exist in the case of a general AP structure on 
the manifold X, see [35] for the proof and construction of the basic contact forms. 

Example 3. Basic example is, of course, the space-time decomposition where 
X = T x B is the product of the classical time axis and the 3-dim space manifold 
B (material or physical) and where T{X) = (d t ) © T(B). In such a case x 1 = t 
and x i ,i = 2,3,4 are spacial coordinates. We take K = T(B) C T(X) to be the 
subbundle of derivatives in spacial directions. 

Example 4. Mathematically trivial but very important physically is the case of 
the Rational Extended Thermodynamics, [21] . where the configurational space is 
extended enough so that the constitutive relations do not depend on the derivatives 
of the basic fields y^. To include this case in our scheme we take K = {0}. Then, 
the bundle J^ n y(ir) — > Y has zero dimensional fiber. 

3.2. Space-time splitting bundle Jg(Tr). Employing the construction of the 
bundle J^(tt) for the product structure T(X) = (d t )®T(B) we define now the par- 
tial 1-jet bundle Jg(Tr). We assume that the generic fiber U of the bundle it : Y — > X 
is the direct product of subspaces of fields y % entering the constitutive relation with- 
out derivatives, with time derivative only, with spacial derivatives only and with 
all derivatives respectively: 

[1, . . . , to] = S U S t U S x U S xt =► U = U x U t x U x x U tx . 

Proposition 2. (1) Bundle 7r : Y — > X is the fiber product of the bundles 
tto : Yq — *- X, 7Tt : Yt — v X, . . . with fibers t/o, Ut, ■ ■ ■ over the base manifold 
X :Y = Y xY t xY x xY xt . 

XXX 

(2) The bundle Jg(7r) is the fiber product of affine bundles 

4(tt) - 0(Y ) x J\ 9t) {Y t ) x J\ dxA) {Y x ) x j\Y tx ). 

(3) Bundle Jgiir) — > y is the factor-bundle of the bundle J 1 (7r). 

(4) Partial Cartan distribution Cs is generated by the 1-forms 

uj" = dy" - z?dt, fie S t , uj^^dy 1 *- z\dx A , \i e 5 X ; w" = dy" - z»dx l , fx e S tx 

in the sense that a section a : X — » Jg(7r) is integrable: a — jp(s) for a 
section s : X — > Y if and only if for all i. a*^ — if restricted to the 
corresponding distribution. 
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Remark 3. Denote by Auts(n) the subgroup of the group Aut(n) consisting of 
the transformations <f> € Auts{n) such that <p preserves the space-time product 
structure and, in addition, </> is the fiber product of automorphisms of the bundles 
7To, 7r t , TT-rTTta;. Construction of the bundle Jg(n) is invariant with respect to the 
transformations <f> 6 Auts{n). 

3.3. Prolongation of vector fields to the partial 1-jet bundles. Below we 
will use the prolongations of the automorphisms (f> of the bundle n whose projection 
<fi to X preserves the AP structure T{X) = K © K' (respectively, an S-structure) 
to the partial 1-jet bundle Jx(n) (respectively to Jg(ir)). Infinitesimal version 
of this procedure is the prolongation of a 7r-projectable vector field £ whose pro- 
jection £ generates a (local) phase flow preserving the subbundle K (respectively, 
S'-structure). In the case of a conventional k-jet bundle J k (n) this procedure is well 
known ( [HI E2 HI] ) ■ In our situation the results are similar. We formulate corre- 
sponding results in the infinitesimal case for the bundle J^in) with an integrable 
AP structure. For more general case of an arbitrary AP-structure we refer to [25] , 

Definition 2. (1) Denote by Xk{tt) the Lie algebra of n-projectable vector 
fields £ in Y such that the projection £ of £ to X preserves the distri- 
bution K C T(M): <pt*K = K for the local flow tf>t o/£. This condition is 
equivalent to the infinitesimal condition C^K C K . 
(2) Denote by Xk®k> (tt) the Lie algebra of n-projectable vector fields £ in Y 
such that the field £ generated by £ in X preserves the AP-structure, |15j 
T(M) = K®K': fa*K = K, 4> t *K' = K' for the local flow <p t oft,. 

Lemma 2. Let the AP-structure T(X) = K © K' is integrable and let (x J , x k ) be 
a (local) integrating chart (i.e. K —< d x j >; K' =< d x u >). Then 

(1) A n-projectable vector field £ = £*(x)d x i + ^(x, y)d y ^ belongs to Xx(n) if 
and only if 

i.e. if the components £, k {x) do not depend on the variables x 3 . 

(2) A n-projectable vector field £ = £?{x)d x i + £ M (x, y)d y v belongs to Xk®k'{^) 
(preserves the almost product structure T{X) = K © K' ) if and only if 

i = e(x)d xi = ¥{x^)d xj + e(x k )d xk , 

Proof. We have [£, d x j] = —(d x j ■ £ jl )d x j 1 — (d x j ■ £ k )d x k. This vector field belongs 
to K if and only if d x j ■ £ fc = for all j and k. The second statement is proved in 
the same way. □ 

Proposition 3. Let the AP-structure T(X) = K@K' is integrable and let (x J ,x k ) 
be a (local) integrating chart . 

(1) For a n-projectable vector field £ = ^{x)d x i + £,^{x, y)d y ^ € Xj<(n) the 
following properties are equivalent 

(a) There exist a vector field t; 1 £ X(J^-(n)) such that 

(i) Local flow of the vector field £ preserves the partial Cartan dis- 
tribution Cok is the Lie field in terminology of [13] J. 

(ii) ttio = £. 

(b) Vector field £ has, in a local integrating chart (x ] \x k ) the form 

e = e{x h )d xj +e{x ki )d xk +e(x\y)d y .. (3.1) 
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In particular the projection £ of the vector field £ in X preserves the 
almost product structure K © K' . 

(c) Vice versa, any Lie vector field on Jk(tt) is the prolongation of a vector 
field £ e XkbK'(Y) of the form (3.1) (see [13] Ch.2). 

(2) In the case where these conditions are fulfilled the vector field is unique 
and has the form 

e = e(x^)d x3 +e(x ki )d xk +e(x^y)d y , + Ue - % (3-2) 

(3) Mapping £ — ► £ is the homomorphism of Lie algebras: [S,,?]] 1 = [C 1 ,?? 1 ] for 
all £,T) £ Xk,K'(k)- 

Proof. Let £ = £ l (x)d x i + ^ fM (x,y)d y ^ + A^9 2 m be a prolongation to the partial 

jet bundle Jx(tt) of the vector field £. Then, condition of the preservation of the 
partial Cartan structure is equivalent to the condition that for all the generators 
W A' — — J2j z jdxi of the contact ideal of exterior forms, CgcOf, = ^„ 1u w k^ f° r 
some functions q^ £ C°°( (ir)). We calculate 

3 

- - X! # w - E 4^). ( 3 - 3 ) 

" 3 

or 

(C* - # + c ''»' + [&> - x 3 - z U h J dxj = E tfw E z 3 dx3 )- 

v 3 

This equality is fulfilled if and only if we have 

' £ M - z^f j = 

\>,x} A j 3i^, x j ~ Vv^j- 

Since neither £ M nor £ J depend on zf first system is equivalent to the requirement 
that both £ M and ^ are independent on x k . Then the second condition determines 
q% and third - — + C^"- Z J — ana - * ne prolongation £ takes the form 

described in the Proposition. □ 

Similar results have place for the bundles J s (tt) defined by a space-time splitting 
S of the fields y*. 

Proposition 4. (1) A vector field £ E <¥s(t) preserves the AP-structureT(X) = 
T(B)® <d t > if and only if f = C(x,t)d x i = £ a (x)8 x a +^(t)d t , 
(2) For any S-admissible ir-projectable vector field £ E Xs(tt) following state- 
ments are equivalent 

(a) There is a vector field £ £ X(Jg(jr)) such that 

(i) Vector field £ G «Y(Jg(7r)) is iriQ-projectable and 7Tio*(^ 1 ) = £, 

(ii) Local flow of the vector field £ preserves the partial Cartan dis- 
tribution Cos at Js( n )- 
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(3) Mapping £ — ► £ is the homomorphism of Lie algebras: [CiV] 1 = [C 1 ! 7 ? 1 ] f or 
all£,r)€X K (Y). 

Detailed structure of vector fields £ preserving the ^-structure and form of their 
prolongations £ x to the bundle Jg(ir) is described in [25], Ch.2, Prop. 11. 

3.3.1. Prolongation of infintesimal automorphisms £ £ ^(tt) to the bundles A^(Y), 

i n+(n+l) / T l / \s 
A 2/l ( J W)- 

Definition 3. (Definition- Proposition, |17j . Def.3.3. for a case of euclidian metric 
G) Let a be a pullback to K^Y of a ir-semibasic form a = a J (x,y)r) v on Y. Let 

(1) Then there exist and is unique a vector field £* Q on Af satisfying to the 
following conditions 

(a) Vector field £* a is tta^y y -protectable and 

7TAJY Y *C a = L 

(b) 

C £ , a e^ = da. 

(2) In an adopted chart the vector field £* Q has the form 
C a =£,+ T Q P d P + T Q p *<%, where 

r » = - P (g + -p^+ (g - o^) = - F ^ G( o-^g + ^ G (,); 

(3. 

ffere £ = , 5 = a- 7 (x, t/)9 X 3 . 

(3) Let a vector field £ 6 ^(X) &e ir-projectable. Then the O-lift £*° o/£ coincide 
with the flow prolongation £ defined above. 

Below we will use the lift of a projectable vector field £ to the bundle n+1 ^ Y\ 

Next result allows to lift £ to the infinitesimal transformation on the bundle A^ n+1 ^Y 
leaving invariant the canonical form p IJl uj tl A n. 

Proposition 5. For any ~k -projectable vector field £ € X(tt) there exists unique pro- 
jectable (to Y) vector field £*( n + 1 ) on the bundle Ag" Y that leaves the canonical 
form p a dy a A r\ invariant. That vector is given by the relation 

r (n+1) = c+e°d p „, e° 

where A G = ln(\G\). 



-p^XCi+g-r )-Pv Ty , 



P"— = -PadlV G {i)-p y — , 

(3.5) 
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Proof. We have, for a vector field of the form £*( n + 2 ) = £ + ^ p<T d Pa 

C^+Div^dy 11 Arj) = 

= (di^ (n +i ) +i i Hr l +i)d)(p f _ l dy tl Arj) = dlp^rj-p^C dy^ Arj^+i^n+^idp^Ady^ Arj) = 

= [fdpn A T] + p„d£" Arj — f dp„ A dy^ A % -p^d? A dy" A m + p X G ,idy» A rj] + 

+ [e"dy^Ar 1 -^dp tt Ar 1 +Cdp tt Ady^Ar li ] = [^dp^Arj+p^dy" Ar 1 -Cdp tt Ady^Ar] i + 

+ P& A dy» A v+Pf.C^Gady'" Ar]} + [^"dy" At]- ^d Pfl A n + f dp„ A dy" A m] = 

= W" + P„C„, + Pn(€,i + C^G,i)]dy^ A r\. (3.6) 

Here we have used the relation dr\i = XciV- 

Equating the obtained expression to zero we get the expression for £ pM as in the 
Proposition. □ 

Combining the last result with the prolongation £*° from the Definition-Proposition 
10 and with the prolongation £*( n+1 ) from the previous Proposition and using fac- 
torization by A^ + '" +1 'y we get the following 

Corollary 1. For any projectable vector field £ G X(ir) there exists unique pro- 
jectable vector field in the space A^ + ^ n+1 ^Y - prolongation of preserving the 
n + (n + 1) form prj + p^dy^ Arji + p v dy v A rj. 

Vector field £ x * projects to Ay* 1 )+(™+ 2 )y f tfe vec tor field f 1 preserving the 
form p^dy^ A T]i + p v dy v A rj mod rj. 

Below we will also need to define prolongations of projectable vector fields £ G 
X p (tt) to the space of the n + (n + l)-forms A™^™ +1 ''(J 1 (7r)) on the 1-jet bundle 

J 1 (7r) — > Y. Next statement, whose proof is presented in the Appendix 1, describes 
these prolongations. Local coordinates in this bundle corresponding to a fibred 
chart (W, x l 7 y^) is defined at the presentation of elements of the fibers: 

prj + p>*dy» Ar/i + An + g* A 77. (3.7) 

Proposition 6. For any projectible vector field £ = £ l c\ + + £f 3 z e G X(tt^) 
on the space Jp(7r), there exist unique vector field £* G A^Aj -1- (J 1 ^))) such that 

(1) Canonical n + l-form Q n+1 = (q^oj^ + ) A n (see ()) is invariant with 
respect to the flow of vector field £* ; £j » Q n+1 = 0. 

(2) This vector field is given by the following expression 

€ = -qffiy* -q v div G (Q)d qv +(-Q&3 -&i, z » -qidiv G (0)d d . (3.8) 

3.4. Partial higher order jet bundles Jp(n). Here we introduce the higher 
order partial jet bundles defined by an AP-structure T(X) = K © K' and, more 
generally, by an S'-structure. In the rest of the paper we will understand 

by Jp (it) , k — 1 , . . . , 00 one of the bundles J k (w) , Jjj- (n) , J| (7r) , J (ir) . One can 
introduce more general notion of higher order (and infinite) partial jet bundle of a 
bundle ir, but this will be done elsewhere. 

Let Z'l be the set of multiindices {I = {h, . . . , i n )\i t G N}. Let T(X) = K © K' 
be an integrable AP-structure on X, and let U, (x 3 , j = 1, . . . x k , k = nj< + 
1, ... n) be a (local) integrating chart in X. Distribution K (respectively K') in this 
chart is generated by vector fields d x j (respectively, by vector fields d x k). Let Nk — 
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{(l,0,...,0),(0,l,...,0)...(0,...,0,l njf ,0,...) C ZD is the set of multiindices 
corresponding to the derivatives in K. Finally let Nk — Nk + Z™ be the set of 
all multinidices larger or equal to the multiindices from Nk with respect to the 
natural order in Z 7 ?. 

Definition 4. Let ttxy ■ Y — > X be a fiber bundle and let K C T(X) be a 
subbundle. 

(1) Let x € U C X, S\,S2 G r(f, 7r)|si(or) = s 2 (x). Sections si,s 2 are called 
K-equivalent of order k at a point x G X : s\ s 2 if d 1 s\{x) = 
d I s 2 (x), VLeN K . 

(2) J K x {tt)- space of classes of ~ K at a point x G X . 

(3) J K (ir) = ^-ixexJfc x( n ) ' the space of partial k-jets of sections ofir. 

(4) Infinite jet bundle ("") is the inverse limit of the bundles J K (tt) under 
the natural projections TTuk-i) '■ J ski 1 *) ~ * 4 _1 ( 7r )- 

Speaking simply, bundle J K (ir) carry information about the derivatives d x j s of 
section s G r(7r) and all derivatives d I d x jS up to the order k. 

In the next proposition we collect some basic properties of bundles J K {ir) that 
will be used below. 

Proposition 7. (1) Bundles ~Kk(k-\) '■ 

J k (tt) J^ _1 (7r) form the inverse se- 
quence of the affine bundles. 
(2) There is a canonical surjection of affine bundles wk ■ J k {^) —* Jxi 71 ) suc ^ 
that the diagram 

J"i*) 4W 



"■fc(fc-i) 



"■fc(fc-i) 



I 



7 (tt) ► J K (tt) 

is commutative. 

(3) Lef T(X) = K ® K' be an integrable AP-structure, {x^ ,x k ) - local chart 
integrating the AP — structure. The 1-forms 

f w" = dy» - J2 jeK z^dxi; i G T^h~, 

[ul = dz^ - J2i zjidx 1 , J e N K , < | J\ < k, fi G T7m 

generate the partial Cartan distribution C K on the bundle J K (n) —> X . 

(4) The "total derivative" operator di — c^+X^i ^ eP ^ eP Zj +i d z » 
is defined on the functions C°°(J k (tt)) and has the same properties as the 
usual total derivatives (commutativity, etc). Here and below we use the 
notation {jJi,I) G P for the set of variables Zj, \I\ > that are present in 
the partial k-jet bundle Jp{^)- The set P if indices is invariant under the 
admissible automorphisms of J k (tt) 

(5) A section q of the bundle Jk( 7T ) * s the (partial) k-jet of a section s : X — > Y 
if and only if q*(ujj)\K = for all 1-forms listed in p. 3. 

(6) Any ■n-projectable vector field £ G X(Y) preserving the AP-structure T(X) = 
K © K' can be uniquely prolonged (by the flow prolongation) to the pro- 
jectable vector field £ k G X(J k (k)) preserving the partial Cartan distribu- 
tion C K . For a vector field £ = Cd x i + i^d y n the k-th order prolongation 
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has the form (see |24j . Thm. 2.36) 

(7) Contact decomposition (2.3) of the forms is valid for the tower of partial 
jet-bundles J^(tt). 

Proof of almost all statements in this Proposition is straightforward, using the 
standard structural properties of the higher order jet bundles, see [12l [24]. In the 
last statement the prolongation is the standard flow prolongation of an admissible 
vector field, see Proposition 3 above, and Lemma before it. 

For a space-time splitting S of the space of fields y l , see Sec. 3. 2, one can similarly 
define the bundles J|(7r), J|°(7r), their partial Cartan structure and establish the 
prolongation properties of 7r-projectable vector fields £ £ X(Y) preserving the S- 
structure of the bundle 7r and the existence of the contact decomposition of exterior 
forms. We will use these results without further references. 

4. Constitutive relations (CR) and their Poincare-Cartan forms 

Here we introduce the constitutive and covering constitutive relations of order 
k as generalized Legendre transformation from the (partial) k-jet bundle Jp(7r) to 
the extended multisymplectic bundles of n + (n + l)-forms on the manifold Y, see 
next commutative diagram 



n+{n+\) 
2 




X n 

Definition 5. (1) A covering constitutive relation (CCR) G of order 

k is morphism of bundles C : Jp{^) — > A% n+ Y over Y: 

C(x\y»,z>t, = (x\y»;p; fj; H M ), p, f£, U, £ C°°{J k p (tt)). 

The Poincare- Cartan form of the CCR C is the form 

©c = C "*(©2 + ©2 +1 ) = PV + Kdy 1 + IW A V- (4-2) 
(2) A constitutive relation (CR) C of order k is a morphism of bundles 
C : J*(tt) -f A£+ (n+1) Y over Y: 

C(x\ y», zf , . . . , z^ Ak ) = (x\y^ F^), F* , II M £ C°°(J p fe (tt)). 
The Poincare- Cartan form of the CR C is the form defined mod r\ 
e c = C*(©2 + ©2 +1 mod V) = F ldy^ A m + Tl^dy" A n mod n. 
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When the order of a CR or CCR is one (k = 1) we will omit the words "of order 

1". 

Any covering constitutive relation defines the corresponding constitutive relation 
by the projection - C = q o C. On the other hand, there are many ways to lift a 
constitutive relation to the CCR. One of them (for k = 1) is to use the connection 
on the bundle 7T, see [14]. In between all possible CCR corresponding to a given 
CR there is one privileged (as we will see in the next section) 

Definition 6. The lifted CCR C of a constitutive relation C (of order k) is 
defined by 

O d = F^u" Aru + Tl^u" A 77 = — ( 2 i ^ + A Vi + H M dl^ A 77. (4.3) 

As the next Lemma shows, the Poincare-Cartan form of the lifted CCR can be 
obtained with the help of the vertical endomorphism S v . Proof of this Lemma is 
straightforward. 

Lemma 3. Let k = 1. For a 1-form \ = F^dz? + n M c?y M S A 1 (J 1 (7r)), one has 

e d = s;(F^+n li dy> i ). 

Remind, that the Poincare-Cartan form Ol of the the Lagrangian form Lrj is 
an example of the Lepage form corresponding to the Lagrangian L (ref). Having 
defined the Poincare-Cartan form of a balance system with the covering constitutive 
relations C it is interesting to see when such a form may be a Lepage form and if 
yes, to which Lagrangian it corresponds. Next result provides the answer to this 
question. 

Lemma 4. The term 0^ of the Poincare-Cartan form (9.4) of a covering consti- 
tutive relation C is the Lepage form in J k (tt) if and only if C is the CCR of order 
one, i.e. that p, F^ € 7rj c *C°°(J 1 (7r)). If this condition is fulfilled, the associated 
Lagrangian of the Lepage form & c is equal to L = p + Y^t^ i)ep z i^- 

Proof. To prove the first statement we write the form O^, as follows 

Q n e = pr, + F^ Ar/i = (p + zfF^r) + F*uf A 

and denote p = p + z^F^. Now we calculate 

Trf +1) *d9£ = 7rf +1) *(# A 77 + dF^ A u" A m + fjdw" A m - F^X G ^ A 77) = 
= (d h p + d v p) A 77 + (d h F* + d v F 1 ^) Aw"A 77, - A dx 3 A m - A 77) = 

= ( P,z^)^+djF^ A^Arn+i ]T F^^A^Am-F^An-F^Xc^An^ 

(ia/),|/|^0 (i/,/),|J|^0 

= -(diF^+F;X G ,i)^Ar]+p, y ^ u An+(p^-Fi)^Ar]+ P^^A^Con. 

In the last expression first four terms present 1-contact part, the last one - 2-contact 
part r u | 7 |> z vLuf) Alu^ A rji. 1-contact part is 7Tq +1 -horizontal if and only 
if conditions 

= 0, |/| > 1, 
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First condition requires that p £ 7rJ c *C 00 ( J 1 (7r)). If this fulfilled, then the second 
condition requires that £ 7rj c *C°°( J 1 (7r)) and therefore the CCR C is of the 
order one. What is left is that the second condition tells that the density/flux part 
components of constitutive relations come from the "Lagrangian" p. Below such 
constitutive relations will be called " semi-Lagrangian" . □ 

4.1. Examples. Present now several examples of different types of constitutive re- 
lations or covering constitutive relations and corresponding Poincare-Cartan forms. 
In all these examples k = 1. 

Example 5. A Lagrange constitutive relation of order 1 is defined by a 
function L £ C°° {J 1 {%)): 

Fit fit 

C L { X \y^z») = {pl = F; = — r ^ = —). 

Balance system for such a constitutive relation (see next section) coincide with 
the conventional system of Euler-Lagrange Equations (2.10) of the first order La- 
grangian Field Theory. 

Remark 4. For k — 1 Euler-Lagrange Equations (2.10) have the form of a balance 
system with canonically defined densities and flux components. For k > 1 Euler- 
Lagrange equations can also be written as a balance system but in this case there 
are different ways to specify flux and even density components due to the presence of 
higher derivatives. This non- uniqueness is similar to the non-uniqueness of a choice 
of Lepage form for higher order Lagrangian Variational Theory, see [4] . That why 
in the Lagrangian case we take k = 1 only. 

Example 6. A semi-Lagrangian CCR Cl,q is defined by a functions L, Q M , /j, = 
l,...,m £ C^iJ 1 ^)): 

BT 

C L ^{x\y», zf) = {p = L~ z?L tZ »,pi = — ; n^Q^x\y^ zf )). 

Example 7. L + £>-system. Let L Lagragian, D (" dissipative potential" ) £ 
C°°(Jp (-71")). Let the time derivatives of the fields j/ M - Zq are present in the partial 
1-jet bundle Jp(7r). The CR Cl,d is defined by its Poincare-Cartan form 

Ql.d = L^dy" Ai]i + (D z » o - L^dy** A r\. 

balance system for such constitutive relation have the form 

SL _ 3D 
Syt 1 dyv 

which is well known in the Continuum Thermodynamics, [18[ I19j . 

Remark 5. If we would like to define directly the analog of Lagrangian or semi- 
Lagrangian relation on a partial 1-jet bundle we would have a situation where an 
absence of a variable from the partial 1-jet bundle leads to the nullity of the 
corresponding flux component. This restricts an application of a conventional semi- 
Lagrangian CR defined on a partial jet bundles. On the other hand, it is possible 
that the CR are partly variational in the sense that part of the flux components 
are defined by a "partial Lagrangian" while other components are to be defined 
independently, see below, Prop. 10. 
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Example 8. Consider a case where only spacial but not time derivatives of fields 
y^ enters the constitutive relations: Jp{^) = Jt(b)( % )' Then it is possible that for 
a function L E C°°(J^ {b) {tt)), = J^r, A = 2,. .. ,n; IT M = §^ but densities 
F® are to be defined independently. Then, 

dL dL 
9 C = F°dy» At]q + Q^dy» Ar) A + g^dy" A V . 

The system of balance equations corresponding to this CR is, of course, first order 
by time derivatives. 

Five fields fluid system (Sec. 2. 7) is an example of such balance system. 

Example 9. Vector-potential CR. (RET case, dual variables) In this case 
Jp(ir) = {•} has a zero dimensional (point) fiber over Y. Let h = h l (x,y)rii be 
a semi-basic (n-I)-form on Y. Define the CR 

dh l 

C h {x\y») = (pl = —-,n» = n iJi {x 1 y)). 

4.2. Form Kc and the Helmholtz coindition. Let C be a constitutive relations 
of first order with the lifted Poincare-Cartan form 0g = + ILjW M A r\. 
Associate with this constitutive relation the following contact (n+l)-form on J 2 (ir). 

K c = (i>f + I^w") A V . (4.4) 

Calculate differential of this form: 

dK c = [(dyuU^Ady^ + idy.Ff-d^n^dy^Adz^Hd^F^dz^Adz^Ari. (4.5) 

Now we use the following form of Poincare Lemma 

Lemma 5. Let fi,i = 1, . . . ,k; gj, j — 1, . . . , s be functions of all variables z l , i — 
1, . . . , k; y J , j = 1, . . . , s such that the form K = fidz 1 + gjdy J is closed: dK = 0. 
Then, locally, K = dL, i.e. fi = d z iL, gj — d y jL for some smooth function 
L = L(z\yi). 

Applying this Lemma to the second form of equality (4.5) we get the following 
version of (local) Helmholtz condition (comp. [?]) for a balance system of order one 
to be Euler-Lagrange system for some Lagrangian L 6 C°°(J 1 (7r)). 

Proposition 8. For a constitutive relations C of order 1 the following properties 
are equivalent: 

(1) Form Kc is closed: dKc = 0, 

(2) C is locally Lagrangian constitutive relation: C — Cl for a (locally defined) 
function L £ C 00 ^ 1 ^)). 

Remark 6. Construction of the form Kc is related to the mapping uj^ f\r\i — > ^ An 
of J 1 (7r) — > J 2 (tt). It would be interesting to construct such a mapping from J k (ir) 
to J fc+1 (7r) for k > 1. 

4.3. Variational sequence and the balance systems. Condition of the Propo- 
sition 8 for a balance system to be Lagrangian points to the possibility to interpret 
and study the balance systems using the variational bicomplex, see [U [24] , or the 
variational sequence, [29] . In this section we realize this possibility. In presenting 
the variational bicomplex we will follow [29 . 

We will be using the augmented variational bicomplex in the form 
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R 



R 



4 



d 4 



0.1 



n n (M) 



-1 



d„ 



1,1 



> f2 n_1 (M) -^=-> — Bq 1 '"- 1 



0,n d„ 



l.n 



0,n ei 



l.n 



d„ 



E p ^ 



ih\ { - 1)Pdh \ { ~ 1}P+ldh [ 



Tp P ,l dv 



E p 



-i.i 



iTip+l,n— 1 
-^0 



ij (_1)Pd "| (-l) p+1 dJ 



d„ 



^0 



+ l,n 



+ l,n 





(4.6) 

In this diagram, E®' q = i7 0,9 ( J°°(7r)) - the bundle of smooth horizontal q- forms 
on J°°(7r), Efr q = C p A^ ' 9 (^°°( 7r )) is the bundle of p-vertical and g-horizontal 
p + (/-forms. Terms of the lowest line represent the fist derived complex 

E p,n = E^/d^EP' 71 - 1 ) = c p [\n°> n /d h {c p f\^ n ~ 1 ). 

The mappings e\ : E p ' n — > E p+1 ' n , ei([a]) = [d v a] participate in the last horizontal 
line of the bicomplex as well as in the second half of the induced variational sequence 
(where Euler-Lagrange mapping £ — e\ o7r, 7r being the quotient projection defining 



- R £ °'° El'° ^ 



n,0,n-l <4 _o,ti £ r-il.n ei E ,2,n 
^0 — y ^0 - * ^1 — *■ ^1 ■ 



(4.7) 

Operator / is the interior iJwZeroperator / : Q, s ' n ( J°°(7r)) — » f2 s '™(J°°(7r)), defined 

by OB) 

I(cj) = A io z uj - d M (ia < to + dfj, lfl2 (id i w) — . . . . (4.8) 

S L s ^i^2 

Interior euler operator is closely related to the Euler-Lagrange operator £ , namely 
(see PP, Ch. 2), for any Lagrangian n-form A = Lr] 

E(X) = I o d v (X). (4.9) 

That is why the next proposition is hardly surprising 
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Proposition 9. Let C be a CR and Kq = HiU) z A i] + F?uj 1 A -q - corresponding K- 
form introduced in (4-4)- Then, the balance system is equivalent to the equation 

fsi e I{K c ) = 0, (4.10) 

for otff e V(tt). 
Proof. 

I(K C ) = J A [n iV - d^Ffrj)} = J A [Ui - d^Ff 1 - F?\ G J V . (4.11) 

Here we have used the form n — ^/jGldx 1 A ... A dx n . 

Applying i^i for any variational vector field £ e r(V(7r)) and taking pullback by 
the mapping J*s we get equation () in the form 

w'(0 ° .^K^ + f?^g,h - n«] o j°° sr? = o. 

Fulfillment of this equality for all (or many enough) variational vector fields £. is 
equivalent to the statement that s is the solution of the balance system □ 

5. Variational form of balance systems 

In this section we present the invariant variational form of balance system and 
the separation of this invariant form into the m separate balance laws by indepen- 
dent variations. In difference to the Lagrangian case, in general, one has to put a 
condition on the variations that can be used for separating equations. In cases of 
semi-Lagrangian, order 1 or RET balance systems that condition is fulfilled for all 
variations as in the conventional Lagrangian theory. In general this situation can 
be remedied by modifying the source term li^dy^ A i] - addition to it some contact 
form (see below, Sec. 5. 2). We prove the main result for the case of a constitutive 
relation of arbitrary finite order, so, for simplicity we will consider that a CR C 
is defined on the infinite partial jet bundle J£°(ir) but is 7r£°-projectable for some 
k < oo. 

5.1. Invariant form of balance systems. We start with an arbitrary covering 
constitutive relation C and change the sign of source term, i.e. we consider 

6^ = PV + F^dy* Arji- Tl^ A V . 

For a section s e Ty (n) , V C X we request the fulfilment of the equation (Invariant 
Balance System, shortly IBS) 

3 k (snUdQ e _) = 0, £ e X{J k p {*)) (IBS) (5.1) 

for a large enough family of variations £ - sufficient for separation of individual 
balance laws (see Def.7 below). 

Thus, we take the n + (n + l)-form of the form (4.2) and apply first the 

Con-differential d and then for a vector field £ = &dj + £ v d y v + J2\i\>i ^i^zf- 
Adopting the summation by repeated indices we recall that only Zj or derivatives 

by these variables with (/i, I) G P are present in the formulas. We introduce the 

notation dr\i — \o,m- It will be convenient to include variables y M into the family 

of variables Zj for |7| = taking z l — y l . 

We will also use the contact splitting (2.3) of the lift of differential of a function 

p to J£° (ir) and similar contact decomposition for the flux components F^ . 
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We have, 

dB d _ = d(pn+F*dyi t Arii)+IL lt dy> t Ari = dpA^+dFiAdy^Arji-F^dy^AXa^n+Il^dy^Ari 

(5.2) 

Now we continue the calculation replacing in the last two terms in (5.2) the dy^ by 
uj^ and using dy v A r]j — u v A r)j + z v k dx k A rjj — uj u A r)j + z^rj 

= [{d lP )dx l - J2 P,z^i] A7 l+ 
(pJ)eP 

+ [(d i Fi)dx i - (Fi) tZ ^}Ady»Ar 1j -X G , i FXAr ] + Ii^Ar ] = 

(^J)eP 

= {d i Fi)dx i Ady^Arij-X^F^An + Il^Ari-i E p^]A V + 

-[ E (Ft), z ^}A(ujiAr lj + z»ri) = 
(p,i)eP 

= (-^F^Av-XciF^Ari+U^Av- E k^+^( F i),^H Av+2Con = 

(pJ)eP 

= [-d tl F?)-X G ^F?+Il l ]u; l A7 1 + F J "LuiA V +[ E (P+^ F l),z^i M + 2Con. 

(pJ)eP 

(5.3) 

Here we have used the equality (d i Fl)dx l Ady' J, Ai] J = -(djFi)dy^ An = -(djFl)uj^A 
rj. Term 2Con in the last formula represents 2-contact form E( M i)ep( F v),z >lU, i] A 
uji A rjj . 

Last formula proves the first statement of the next 

Theorem 2. Let C be a CCR defined in a domain of the partial k-jet bundle Jp(jr), 
k ^ 1. Then, 

(1) We have the following contact decomposition 

de c = [-dify - Xg^I + n M K a v + a v - 

E (P+ E ^),*? W J A77 + 2CW (5.4) 

Term in the brackets is the vertical differential d v (p + ZjFj>) of the func- 
tion p + ZjF-j,, see J5J[T3]. Internal sum is taken over all indices such 
that z» € JpV) 

is in the Tjet projection of j£ (tt), while the outside sum is 
taken over all (p,I) with Zj in Jp(ir). 

(2) Let £ G X(Jp(n)) be any vector field. Then, 

itdB d _ = -^ g (0v - " h 5 +1 (t k+1 )v + Con, (5.5) 
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for an arbitrary prolongation £ fc+1 to the vector field £ to Jp +1 (7r), where 

= E(^ e p^>;(^ +1 ) - [E( ft i) e p(p+E W)e p^),^?^ +1 ) 

(5-6) 

-ffere £ 2 is the projection of £ to J%(ir) for k ^ and f/ie prolongation to this 
subbundle for k = 1. 

(3) For £Ae Zi/fed CCi? C of a constitutive relation C and an arbitrary vector 
field (tt)), 

i f d0 e _ = -4(£)ry - E F^(e +1 )v + Con (5.7) 

(4) Decompositions (5.4,5.7) are covariant with respect to a change of admissi- 
ble coordinate system (x l ,y^) — > (x n = h l (x 3 ) , y' M = h tJ, (x J ,y u ) preserving 
the "partial" structure (both in K (& K 1 and S cases). 

Proof. Previous discussion proves the validity of the first decomposition. Second 
and third statements follow from the first one. 

Last statement follows from the tensorial behavior of both components of de- 
compositions with regard of such coordinate changes. □ 

Remark 7. Quantities = w M (£) = £ M — zf£ l form the characteristic of the 
vector field £ = t^di + £ M <9 yf i in the sense of [24j . Ch.2, or the generating section in 
terms of |13j . 

Remark 8. Equality (5.7) contains the jet variables of the second order zfj. Yet, 
only u>? with (i,fi) G P are present in the formula (5.7)! For instance in the 
RET case all these terms are absent from (5.7). 

It is seen from the formula (5.7) that the lifted covering constitutive relation is 
privileged in the sense that the formula for dQ^ simplifies essentially for lifted 

CCR CL. 

In order that the equation resulting from taking the pullback by j k (s) in (5.4) 
would not depend on the variables not in J^tt) we require that all the 
coefficients of these variables would be zero. This leads to the condition one has to 
put to the allowed variations £ (see Def.7 below). 

Consider now two cases where no restrictions to the variations of Poincare-Cartan 
form appears. 

Proposition 10. Let C be a CCR of order k. Then. 



e f h 



E *i J ■■■'■■) 

(pj)eP Oj)eP 



(5.8) 



if and only if for some L G C°°(Jp (it)), F* = d z ^L, {fi,i) G P. i.e. CR C is 
(locally) semi-P-Lagrangian (lagrangian by variables z^\(i,fi) G P ) and 
does not depend on the variables Zj , \I\ > 1. In this case, p = £— y\ z^d z ^L+ 
l(x,y) with an arbitrary l(x,y) G C°°(Y). 



22 



SERGE PRESTON 



Proof. Combining terms with |/| = 1 in the sum with the first term in (5.8) and 
using linear independence of basic contact forms we split condition (5.8) into the 
following two groups of conditions 

fl. Fi - %(p + £ M)eP z£F*) = 0, V (j, v) e P, 
\2- (P + E^jeP^),,? = 0, V ( M , J) G P, |7| > 1. 
Rewrite the equalities of the first group in the form 

d z »p = -z£d z »F*. (5.9) 

To sec that => holds we notice that provided left statement is true, the right sides 
of these equalities satisfy to the mixed derivative test 

d z »(-z?d z »Fi) = d z »(-z?d z »Fi), 

or d z »F 1 ^ = d z ^F^. It follows from this equality valid for all couples of indices 
(v,k),(n,l) G P that there exists a function L G C°°(Jp(7r)) such that Fl = 
d z vL,\/(v,j) G P. Substituting this to the first family of equalities we see that the 
function p + £^ k ^ eP z%L z g — L does not depend on z?, (fx, i) G P. Therefore, 

P = L- 4 L ,z-+l(x,y,zr, (n,l)£P). 

(n,k)eP 

Moving sum in the right side to the left and using second group of equations (5.9) 
we see that the function L + 1 does not depend on the variables z", (y, I) € P \ P. 
Using this function instead of L we get the conclusion of Proposition. 

To prove the opposite - reverse the arguments. □ 

Notice that the conclusion of this proposition put no restrictions on the compo- 
nents Ft*, (fi, i) £ P. 

In the case of a semi-Lagrangian CR there are no restriction to the variations in 
the IBS equation (5.1) provided one uses the proper lift of the constitutive relation 
to the CCR (see previous section, Example 2). More specifically, 

Theorem 3. Let Cl.ii be a semi-P-Lagrangian (by flux components in P) consti- 
tutive relation of order 1: F'^ = L z ^, (n,i) G P where L G C°°(Jp(ir)). Let Clji 
be a CCR covering C with p = L — J2(pi)eP z i^z^P- Then the following statements 
are equivalent 

(1) For a section s G T(tt) and for all £ G X(Jp(n)) 

j 1 p ( s yi i dG dLii _=o. 

(2) Section s is the solution of the system of balance equations 

i\(n,i)eP MtP 

(5.10) 

Remark 9. In a case of a Lagrangian Field Theory of higher order (k > 1) the 
Euler-Lagrange system of equations does not have unique representation as a bal- 
ance system. Perhaps this is the reason why the last results are limited to the 
Lagrangian system of the first order. Probably introduction of more general form 
of balance systems allows to include higher order Lagrangian systems in this scheme 
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in such a way that the arbitrary vector fields £ in the fc-jet bundles are admissible 
as variations. 

Similarly, there are no restriction to the variation in the RET case where CR is 
the section of the bundle 7r™+(™ +1 ) : A^ n+1) Y (see [25]). 

Theorem 4. Let C be a constitutive relation of the RET type and C - corresponding 
lifted CCR. Then, for a given section s £ r(7r) 

s*(i^dQ^ ) = V £ € X(Y) Section s is the solution of -fa. 

In general case we have to introduce the class of admissible variations. 

Definition 7. A vector field £ s X(U), U C Jp{^) is called C -admissible <E 
*c(U)) if 

(1) In the case k = 1, for some (—any) prolongation of £ to the vector field 

ezx(j*(*£(u)) 

(2) For k > 1, for some (=any) prolongation of £ to i/ie vector field g 
^(^ +1 (^( fc 1 + i)i(^)) 

<4 +1 (0 = E (tf-* fc 4)*j = °- 

(M,i)eP 

Remark 10. Condition of C-admissibility is much more restrictive in the case 
k = 1 then in the case k > 1 since in the last case may depend on z? k . 

As the next result shows, locally any CR is separable. 

Lemma 6. Let W C Y be the domain of a fiber chart (x l ,y^). Any CR C is 
separable in W k = 7r^" (W), more specifically, Vy £ W, z e W k , TTko(z) — y, the 
mapping Xc{W k ) 9 £ — > G * s ^ e epimorphism. 

Proof. Consider vertical vector fields that are constant along the fiber U fl W. □ 

As a result, in general case we have the following result 

Theorem 5. For any CR C and a domain U C Jp{^) the following statements for 
a section s S T(jv)(U), U C X are equivalent: 

(1) C is separable in U and for all £ £ T(U,Vc(Jp(^))> 
i fe ( S )*(i s de e J = 0. (755) 

(2) Section s £ r(f , 7r) is i/ie solution of the system 
(F;oj k (s))^+F^oj k (s)X G ^U^j k p (s)), M =l,...,m. (*) 



24 



SERGE PRESTON 



5.2. Contact source correction of the lifted PC-form 8^ and the balance 
equations. In this subsection we show that one may modify the Poincare-Cartan 
form 6p by adding a contact term to the source in such a way that applying the 
procedure of variation described above to the PC-form modified in this way one 
can remove the restriction s ont the variational vector fields. 

Let A2 + ^™ +1 '(Jp(7r)) — > Jp{it) be the bundle of n + (n + l)-forms on the partial 
1-jet bundle Jp(it) annulated by substitution of any two ir 1 - vertical vector fields. 
An element of a fiber of this bundle have the form 

PV + P^dy* Arji + p^dy» Arj + k^dz^ A Tjj + q^dz? A rj. (5.11) 

Introduce on this manifold the following n + (n + l)-form: 

Q = prj+p^dy^ A rj, +p ll dy f * Arj + q^dzf Arj = pij+p^dy 11 A rj t +p^ A rj + q^uf A rj 

(5.12) 

which differs from the canonical form G>2 + 2 l+1 by the last term only. It is easy 
to see that this form behaves covariantly under the change of fibred coordinates. 

Let C : Jp(vr) — ► h^ + ^ n+1 ^ (Y) be a covering constitutive relation of the form 
C(x, y, z) = prj + F^dyV Afjj- Ii^dy^ A rj. Mapping C can be lifted to the mapping 

C 1 : J k p (7r) - A^ +( " +1) (JpW) : (x,y, z) - pr] + F^dy* A Vi - K c , 

where the "source" form Kq = H^u)* 1 Arj + F^w? Arj was introduced in the previous 
section. 

Now we will modify the Poincare-Cartan form 6^ on the (partial) jet bundle 
Jp (ir) for fc > 1 and on the partial 2-jet bundle Jp{ii) for fc = 1 by adding an extra 
source term to get the modified covering Poincare-Cartan form 0^ by the formula 

Q d = C 1 *Q = prj + F^dy^ Arji — Arj — F^ A rj. (5.13) 

Applying to the form 9^ arguments leading to (5.4-5.7) we get, for a lift £ 2 to 
the bundle J%(n) of a vector field £ G X(J^(n j) for fc = 1 and for a lift £_ k+1 G 
X{J k+1 {i:)) of the vector field £ k on J k (n) for fc > 1 

; £2 Je£ = {-u/(0[^+A Gli F^ 

+ FX{e)rj - i>f A i S Tj + Con = 
= {-^(0[^+A Gii F;-n^-^(^ 

(5.14) 

since the second to the last last form in the previous expression is contact. For the 
lifted covering constitutive relation C (see Def.6) this expression simplifies: 

i e dQ d = -^{^[diFl + X g ^F; - n„]jy + Con. (5.15) 
Thus, we have proves the following 

Theorem 6. For a given constitutive relation C of order fc and a section s G L(7r) 
of the bundle ir the following statements are equivalent: 

(1) For any vector field £ G X(Jp(n)) and its arbitrary prolongation to the 
n k -projectable vector field £ fc+1 on J k+1 (Tr), 

j k+1 *(s)i ik+ idQ e _ = 0, (5.16) 
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(2) Section s is the solution of the system (F^ o j k {sj) = n AI oj' c + 1 (s), fi = 
l,...,m. 

5.3. Dual form of a balance system. Reversing in the IBS equation (5.1) the 
order of contraction with a vector field £ fe and acting by the Iglesias differential 
(but keeping the sign of the source term) we get the equation 

/ +1 (s)*<%0 c = 0. (5.17) 

Calculations similar to one performed above for the equation (5.1) allows to prove 
the following 

Proposition 11. Let C be a CCR and let C = prj + Ffdy 1 A + U t dy l A rj 
be the corresponding Poincare-Cartan form. Then, for any projectable vector field 

<%© c = [u*(0[diFi + f;a g , 4 - ry + + x G , t ) ■ ( P + ]T z»Fp+ 

+ (p+ ]T z^FpdiC + F^^rj. (5.18) 

For the lifted CCR C, where p = — J2(^ i)eP z i^> ^ e P rev i° us formula takes the 
form 

<%6 C = [^{0\diFl + F;A G , 4 - LLJ + F^(0] V (5.19) 

Remark 11. Notice that the condition F^diU^ (£) = is equivalent, for the vertical 
vector fields £ G V(n) to the admissibility condition of Dcf.7. 

Corollary 2. Let C be a constitutive relation of order k ^ 1 and let C be the 

corresponding lifted CCR. Then for a section s e T(U,n), U C X defined in an 
open subset U C X such that C is separable in 7r _1 ([/) ; the following conditions 
are equivalent 

(1) For all C-admissible vector fields £ € 7r _1 (L/) 

j k+1 (s)*di^e d = 0. 

(2) s is the solution of the balance system 

(Fi l oj 1 (s)),.=U^oj 1 (s), M =l,...,m. 

Remark 12. Expression in the left side of equation (5.17) applied to the 0^ 
represents the second, boundary term in the Cartan formula for the Lie derivative 

£ ? C = i c deg + di^. 

It is interesting that one can use both these terms to get the balance system 

6. Geometrical transformations and the symmetry groups of 
constitutive relations. 

Let C be a CCR with the domain Jp(n), let (f> s Aut p (Tr) be an automorphism 
of the bundle n preserving the partial structure. Denote by <f> k the flow lift of <fi to 
to the manifold J k (ir) and by (f) 1 * - its canonical flow lift to the space 
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Define the action of <f> on the CCR C by the "conjugation" 

= 1 * oCo<j) k - 1 . (6.1) 

Similarly one defines the action of <f> on a CR C, using the projection cf) 1 * of 1 * to 
A «+("+i)y instea d of see [HQ7]. Immediate check shows that = (C^)f 

Since the transformation 1 * preserves the canonical multisymplectic form on 
A.J+( n + 1 )y it is easy to see that for the CCR one has 

G 6 , = 4> k ~ u Q d . (6.2) 

6.1. Geometrical symmetries of CR and CCR. 

Definition 8. (1) An automorphism <fi G Aut p (ir) is called a geometrical 
symmetry transformation of a CR C (resp. of a CCR C) if 

=C {respectively, & = C). 
(2) Let £ G X p (tt) be a ix -protectable vector field. £ is called a geometrical 
infinitesimal symmetry of the CR C (resp. of a CCR C) if (see |12j . or 
below for the definition of a Lie derivative along the mapping) 

L,gk£i t -,C = 0, 

(respectively, L(tk£i*\C = 0). 

In the next proposition we collected some simple properties of geometrical sym- 
metries and infinitesimal symmetries of constitutive relations and covering consti- 
tutive relations. These properties follows directly from definitions and the fact that 
the transformation (j) 1 * preserves the canonical multisymplectic form on K^^^^Y . 
In this proposition we put 

W p = J» x K^ n+l \Y),W p = J» x A- +(n+1) (Y) 

Proposition 12. For a CR C (respectively for a CCR C) 

(1) Geometrical symmetries <f> of C (resp. C) form the subgroup Sym(C) C 
Autp(n). 

(2) Infinitesimal symmetries of C (resp. C ) form the Lie algebra q(C) C X p (k) 
with the bracket of vector fields in Y as the Lie algebra operation. 

(3) A vector field X G X(W P ) (resp. X G X(W p )) is the generator of the 1- 
parametrical group of generalized symmetries of C ( resp. C) if and only if 
it is tangent to the graph Tc (resp. Tg). 

(4) A vector field £ G X p (tt) is an infinitesimal symmetry of C if (and only if) 
the (local) phase flow diffeomorphisms $f — ipt x <f>-t of W p (resp. ofW p ) 
defined by the prolongation of £ maps Tc into itself satisfies to the relation 

$«(r c ) = r c resp. ,*|(rg-) = r f . 

i.e. if the (local) phase flow <f>t of vector field £ is formed by the geometrical 
symmetry transformations of C (resp. of C ). 

(5) Let £ be an infinitesimal symmetry of a CR C (respectively of a covering 
constitutive relation C ). Then, for the local phase flow 4>t of £ one has 

e c * t = $ - 1 *e c = e c ^ c e @ c = o, 

(respectively, C^kO^ = 0). 
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In the next statement we collect the basic properties of action of automorphisms 
4> G Aut p (ir), including geometrical symmetries of constitutive relations, on the 
balance systems Be- For a domain U C X denote by Solc(U) the space of solutions 
s : U — > Y of the balance system it. correspondence U — > Sol(U) determine the 
sub-sheaf on the manifold X. Denote by Sole the sheaf corresponding to the 
introduced subsheaf. 

Theorem 7. (1) Mapping £ — ► maps the space of C-admissible vector 
fields Xc(W) at a domain W onto the space of C^- admissible vector fields 
X c <t,((p(W)) in the domain 4>{W) for any open subset W C This 
mapping defines the isomorphism of (pre) -sheaves of admissible vector fields 
Xc ^ Xc* ■ 

(2) Mapping s — > s^ = ^oso^i" 1 , on the sections of the configurational bundle 
7r : Y — > X maps the space of solutions SolciW) in the domain W C X 
onto the space of solutions Solc<t>(<fi(W)) of the balance system Bq4> in the 
domain <f>(W) c X: 

Sol c {W)^Sol c ^{W)). 

(3) Let (f> G Aut p (ir) be a symmetry of the CR C. Then the mapping s — > s^ = 
maps the sheaf Sole of solutions of the balance system it into itself. 

Proof of this Theorem is based on the following statement 

Lemma 7. Let C be a CCR with the domain Jp{ir) and the Poincare-Cartan form 
0(7 = PV + Ftdy 1 A rjfj, + Hidy 1 A r\ and let <j> G Aut p (ir k ) be an automorphism of 
the double bundle Jp{it) — » Y — > X . Then 

g* Q 6 ] =* pr) Ffdy* A ^ U t df A n, (6.3) 

where 

>=b4+(^° 4W, x »J{4r%\ • detJtf), 

< *F? = detJ$)J(^y v (F» o 4)<tP yt , (6.4) 
^U i = detJ(4>)(U j o^)<fP yi . 

Proof, (of Lemma). We notice that = {detJ(4>))rj where detJ((f>) is the Jacobian 
of the (local) diffcomorphism <f> defined by the volume form rj. On the other hand 

4>*V» = 4>*id x „V = i $z 1 d x ^* T l = detJ {4>) J (4>~ 1 T f Jiu 

since $Z x d xV . = J^ -1 )^*. 
Altogether we have 

4i*Qq = 0* \prj + Ffdy 1 Ai]^ + H i dy i A rj\ = (p o • 

+ (if o fad^x, y) A + IT o ^d(^) A 0*7? = 

= p o^-detJ{4>)n + (if o 4>)(4>i xa dx a + <f>\ yj dyi) A [de*J(0) J(^ 1 )^^] + 

+ IL. o 4>(<P x „dx° + <t>i yj dyi) A detJ($) V = [detJ(4>)J(4>-\(F? o 0)^ yj ]dy j A r/,+ 

+ [detJ^^o^)^,]^ A 77 + [po^-detJ^) + (if o ^[^detJ^) J(^XW 

(6.5) 

Splitting the terms we get the result stated in Lemma. □ 
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Proof. Proof of Theorem 7. Transformation <p, generating the transformation 
of a constitutive relations C — > C^, transforms the (l,n)-component Ff of C into 
the corresponding component Ff M = detJ((f>~ 1 )J((j))^;(Fj o ((f) 1 )^ 1 )^] 3 of the 
CR C*, see last Lemma. On the other hand the mapping 0* maps a vector field 
I = f^sf + f <9 yI + + . . . G <Y(J^(tt)) into the vector field with the local 

2^-component having the form 

see Appendix II in [25 . As a result for the value of the pairing we have 

Since detJ^^ 1 ) > 0, expressions for contraction of flux components and compo- 
nents of lifted vector field in left and right sides vanish simultaneously. This proved 
the first statement in Theorem 6. 

To prove the second statement notice that 

= (£(«))* °vj -v.e^ 1 * e c^ = C?»)* = °kfee_- 

(6.6) 

Since the mapping £ — > 0*£ is, by the first statement the isomorphism of the 
pre-sheaves of admissible vector fields, expression in the right side is zero for all 
£ G Xc(W) if and only if the expression in the left side is zero for all vector fields 
£ G Xc*(<t>(W)) and, therefore s G SoZcCW) if and only if G SW c *(0(W)). 

For the proof of the third statement we use the symmetry condition in the 
form presented at the end of Definition 8 and notice that in this case = C 
so the mapping s — > maps the space Solc(W) isomorphically onto the space 
Sol c ((f>(W)). □ 

6.2. Infinitesimal symmetries of the constitutive relations. Here we con- 
sider the infinitesimal transformations of a CCR C and these of a CR C of the 
first order. We start with an arbitrary covering constitutive relation C : z = 
(x^y\z^^(p(z),F?(z),Tl k (z)). 
Let 

i = e(x)d x » + C(x, y)d yi + &(x, y, z)d^ = £ + &(x, y, z)d z% 

be an arbitrary infinitesimal automorphism of the double bundle Jp{^) — > Y — > X. 

Here we introduced notation £ for the projection of £ to Y . 

In particular, let £ = ^ tJ '(x)d x ti y)d y i G X p (ir) be an infinitesimal automor- 

phism (vector field) of the bundle n, i.e. a projectable vector field in Y satisfying to 
the conditions of Sec. 3. 3 for lifting to the partial 1-jet bundle Jp(7r). Let, as before, 

be its prolongation to the projectable contact vector field in Jp(n). Then, as an 
example of vector fields £ one may consider the flow lift (see (3.2)) 

e 1 = e(x)d x . + C(x, y)d yi + (d^ ~ zi^j 9 4 , 
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z], that are 



of £ to Jp(vr). Here summation in the last term is taken over the 
present in the partial 1-jet bundle. In the RET case we do not need to introduce 
any prolongation. 

Projection £ of a vector field £ defines the prolongation £ x * to the projectable 
vector field in Z* 



j^{n+i)+{n+2)y p reserv j n g canonical multisymplectic form(s) 



(See Sec.7) and, the prolongation £ x * to the bundle A* n+1)+(n+2) (J 1 (7r)): 



di 3 



-ti^Ld -<ljdiv G (Z))d q ». (6.7) 



J qk 



We have used here the equality £ l 3 - =0 valid for automorphisms of the double 
bundle (ir) — > Y — > X. It is easy to see that vector field £ 1 * is projectable to 
•^■i/i" 1 +2 Jp(ir). This projection - £** has, in fibred chart, expression (6.7) with 
the third term omitted. 

Let now <\>t be a local flow in JHtt) of the vector field £ and tpt be a local flow 



A 



(n+l) + (n+2) 



(J 1 (7r)) of the vector field £ x * (respectively ipt be a local flow in 



A^-r^-r^y Qf the vector field 



(n+l) + (n+2) 1 
2/1 

Dcrivating by t at t = the expression for the transformed mapping tptoC o 0_ t 
(respectively for ipt ° C o <f>-t) we get the generalized Lie derivative of mapping C 
with respect to the vector fields (£, £ 1 *) (see [12], Chapter 11) - the vector field over 



the mapping C : Jp(7r) — > A 



(n+l) + (n+2) 



y (respectively C : J^(tt) -> Z for CR C): 



%^)C = r°c-a(£). 



(6.8) 



In local fibred coordinates we have for C^^^C the following expression 



£.p + ^ G (£).p + Ff^ r 



£ • n fc + rfw G (£)n fe + + Fft 



3 Qyk 



n 



d, 



(6.9 



for a CR C the term with d p is absent from the expression of L,| ^i.-jC 

Condition that the generalized Lie bracket (6.9) is zero for a CCR C and a vector 
field £ S AT(y), has the form of a system of differential equation of the first order 
for the components of the covering constitutive relation C: 



£-p + div G (Z)p + F, 



tide _ 



o. 



£ • F£ + <Hvg®F£ + - = 0, V M , k, 

£ • n fe + div G (0n k + n,- J£ + F?zl yh = o, vfc = 1, 

£ • F£ + Fjtdiv G (0 + FT?. = 0, Vfc, /x. 



(6.10) 
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Vector field £ = + £ l <9j + Qd z i in these equations acts on the components 
of the vector functions in the space Jp(n). 

If we subtract forth equation of this system from the second one we get 

F ^- F ^- mt ^ = °- (6 - 11} 
This equation may replace the forth equation in the previous system. Yet, for 
£ = £ , £ G X W (Y) this equation is identically fulfilled. To sec this we remind 
(see ()) that in this case £j = d v £ l — z^^ v = ^ v + z l v Z? yi — z^^ and, as a result, 
= 5 kK? v i ~ = S ^y* ~ S U> Substituting this in the equation () we 

will see that it is fulfilled identically. 

Thus, for the case when £ = £ x for a vector field £ £ X p (ir) the use of modified 
CCR instead of the usual one does not add forth condition to the first three in the 
system (30.6) although it adds an extra term to the third equation. 

We can rewrite this system of equations as the system of conditions for a vector 
field £ £ X p (ir) to be an infinitesimal symmetry of the CCR C: 

n k ■ iKzt?» + + F kA + i F k^v G & - + F i:^ v ] = o, 

ni k : [u k ^ + u^ k + n^e] + [<Kv G ®u k + u k ^e] + f?z; yk = o. 

(6.12) 

In the case of a modified CCR (see Sec. 5. 2) last condition contains both flux and 
source terms while for the original CCR last equation contains only source. Similar 
conditions an infinitesimal symmetry of a CR C are obtained from these by removing 
the first subsystem of equations. 

Consider now the form these conditions takes in two special cases of vector fields: 

Example 10. 7r-vertical vector fields. 

For a vertical vector field £ = C^y* £ "K 71 ") (infinitesimal gauge transforma- 
tion) the system (12.18) takes the following form 

/: p. d,,C • • £< I). 

H k ■ { F U d ° + F i^) & + F tv^ = 0) fj, = 0, . . . ,n;k = 1, . . . ,m, 

I": ■ (n fc;4 rf CT + n^) + - + Fte,, yk = o, k = i, . . . , m. 

(6.13) 

Example 11. ^-lifted vector fields, ^-trivial. 

Consider now mathematically simple but practically important case of a trivial 
bundle n : Y = X x U —> X, fix a fibred chart (W; x M ,y l ) and consider the zero 
connection vq in this bundle defined in the domain W C Y of the chart above by 
= 0. The i/ -horizontal lift of a basic vector field 9 M is — with zero vertical 
components. Its flow prolongation to Jp(n) has the same form: <9* =8^. 

Conditions (6.12) for the vector field <9 M to define geometrical symmetry of a 
CCR C takes the form 

I ■ P,^ + Ag,mP = 0, 

HI ■■ + Xc^Fj; = 0, (6.14) 
[III k : n fc ,„ + X G ,^k = o. 
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since g?«vg(<9 m ) = Ag,^- These conditions require that the densities of fluxes, source 
components and p are invariant under the translation along x^-axis. In a case where 
|G|-const these conditions reduces to the simple independence of all components of 
CCR C on the variable 

Next result establish relation between the infinitesimal symmetries of a CCR 
and those of the corresponding CR. 

Proposition 13. (1) Any infinitesimal symmetry £ (respectively geometrical 

symmetry) of a CCR C generate (by the projection A < £ l+1 ^ +< " n+2 ^ — > Ayi 1 " >+> ' n+2 ' > 
the infinitesimal symmetry (respectively geometrical symmetry) of C. 
(2) Let C be the lifted CCR of a CR C. Then the previous mapping defines 
the bijection between infinitesimal symmetries (generalized symmetries I;) 
of C and C . 

Proof. First statement follows from the fact that second and third subsystems of 
the conditions (6.12) do not depend on p. 

Second statement follows from the fact that for the lifted CCR, where p = 
~ z% jJ ? i \ multiplying equation 11% by z^ and sum by k,/j, we get exactly equation 
I. □ 

6.3. Homogeneous constitutive relations. If the state space of a theory con- 
tains enough background or dynamical fields to make the constitutive relation C 
free from the explicit dependence of C on (t, x) e X (general relativity, theory of 
uniform materials and RET (see [3T]) are three examples), then the correspond- 
ing balance system simplifies and while studying it one does not need to introduce 
assumptions on the character of the dependence of the balance system on a space- 
time point. Definition given below is an invariant way to distinguish a class of such 
CR. 

Any local chart x^ in X defined the local (translational) action of R n in X associ- 
ating with the basic vectors the vector fields d x ^ . Vice versa, any n-dimensional 
commutative subalgebra t) of the Lie algebra of vector fields X(U), U being an 
open connected subset of X, defines the locally transitive action of R n in U and, 
therefore, a local chart in a neighborhood of any point in U. 

Proposition 14. - Definition. Let v be a connection in the bundle n. We will call 
a constitutive relation C v -homogeneous if the following equivalent properties of 
the constitutive relation C are valid: 

(1) For any point z £ Jp(n) there exists a local chart in a neighborhood U x , x = 
7r (z) such that the (n+l)-Poincare-Cartan form flc of the CR C is invari- 
ant under the local flows <j>f of the lifts of ^-horizontal vector fields 
£, £ S t} in the neighborhood of y = ttiq(z): 

jO^i^c = mod Con. 

(2) For all £ G f), the v -horizontal lift £ is the infinitesimal symmetry of the 
constitutive mapping C in sense of Definition 8. 

(3) The graph Tc C Jp x Z of mapping C is invariant under the flow generated 
by (flow) lifts of v -horizontal vector fields 

Proof. Equivalence of statements of this proposition follows directly from the Propo- 
sition 11. □ 
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Similar definition and results formulated in the last Proposition can be formu- 
lated and proved for a CCR C. 

Remark 13. In a case where connection v is flat, the association £ — > S; 1 defines 
the Lie algebra homomorphism fj — * Aut^ 1 ) C X{J^). 

Using first the z/-horizontal lift of vector fields from X to Y and then the flow 
lift to Jp(7r) (Sec. 3) one can get the ^-homogeneity conditions of a CCR C, those 
of lifted CCR C or those of the CR C. 

7. Noether Theorem. 

In this section we present the (first) Noether theorem for balance systems corre- 
sponding to the action of a Lie group G C Sym(C) C Aut p (Tr) of the geometrical 
symmetries of a CR C. We denote by g - the Lie algebra of the group G, g* - its 
dual space, £ — > £r € X(Y) - action mapping on the manifold Y and £ k be its lift 
to the partial k-jet bundle J£(ir), see Sec. 3.4. 

7.1. Momentum Mapping and general Noether Theorem. Condition that 
the group G is the symmetry group of the CR C (and its lifted CCR C) has the 
infinitesimal form (see previous section) 

%(e£ + e^ +1 ) = o. (7.i) 

This condition splits into two corresponding to the order of the forms: 



J >c C fc eg, = (i e d + di e )e 7 t = o, 

[C^Q^ +1 = (t e *d + di e )Q^ +1 = 0. 
For the first condition we have, using formula (5.7) of Theorem 2; 



(7.2) 



c%9g = -i^dQ^ = -i ch [dQ e _ - = 

= <4(0+ E H^{i k+l )v A 77) + Can = 

= W "K)[d ! ^ + A G , ! f;-n (J ]+ F^(e +1 )v-^eri + Con. (7.3) 

Let now s : X — > Y~ be a solution of the balance system 23c ■ Taking pullback of the 
last equality with respect to the jp +1 s we see that the first and last terms in the 
right side vanishes and we get 



3 k {s)*di e Ql = dj k {s)*i e Q n G = -(H^e)oj k (s)r,+ 



E ^K(£ fc+1 ) 



oj k +\s)r). 



(7.4) 

Now we recall definition of the multimomentum mapping, see [6j [14] 



Definition 9. A multimomentum mapping J c : Jp(ir) —> A n (X) <E> g* of a CCR 

C is defined as 
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Lemma 8. 

J 6 {*)(&) = "[(P + + Fy(0\Vi + DCon (7.6) 

where DCon is the sheaf of contact forms whose differential is also contact. 
Proof. We have 

i^ Y e n d (z) = i^iiz) - p?m + Fi^m - Fftd y » a mj = 

= (PC + CF^H - F^uf A ruj - (Fftz^Vi - Fftz? Vj ) = 
= Ip€ + - Fl?z* + FftzflrH - Fftuf A Vij = 

= [(P + z?F>)C + F^(0}Vi ~ Fi?uf A Vij . 

Now we notice that not just the last term in this equality is contact but its differ- 
ential is contact as well, namely 

d[F^ j u» A ruj] = Con - F^ j dz% A dx k A %j = Con - F^ j dz% A (Sfa - S^n-j) = 
= Con - F^e « + <M) A {6fa - 5 k m ) = Con F^ S\ - 5ty)r, = 

= Con-F^{z%-z%) = Con. 

□ 

Using the formula (7.4) we obtain the proof of the first statement of the next 
Theorem. Proof of other statements of Theorem is straightforward. 

Theorem 8. Noether Theorem (general). Let C be a constitutive relation 
defined on a partial k-jet bundle Jp(n), let C be an arbitrary covering CCR (with an 
arbitrary p) of C and C - lifted CCR of C. Let a Lie group G C Sym(C) C Aut(n) 
be a geometrical symmetry group of the flux part 0^ of CR C. Then, 

(1) For all £ e g and for all solutions s e Sol(X, Be) of the balance system if, 
d[J d (j k (s)(x))(&)\ = -dj k *(s)(x)[((P + ^Fp^ + F^)) Vi ] = 



= (IU£)oj*(*) 



*>i(# +1 ) 



°f + \s). (7.7) 



(2) If C — C, then the term (p + ZjF^)^ Y i- n the left side vanishes. 

(3) If Qy C Xc{Y) on Y, last term in the right side drops out. 

(4) The same formula holds for a Lie group of projectable Cartan symmetries 
G C Aut p (wko) of the space Jp{^) and its Lie algebra qy C X(J k (tt))- 

(5) If Qy C Xc(Y) and if the system Be is the conservative laws system (i.e. 
if Hp = 0, a = 1, ... ,m), then V£ e g and for all s G Sol(C) the Noether 
conservation law holds: 

d[j k *(s)(x) F y Y {i) Vl ] = ((F^ Y (0)oj k »(*));* = 0. 
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7.2. Semi-Lagrangian and RET cases. In the special cases of semi-P-Lagrangian 
and RET constitutive relations the formulation of Noether Theorem is conveniently 
simplifies. 

Theorem 9. Noether Theorem for Cl.ii, k = 1. LetCL.n be a semi-P-Lagrangian 
CCR with 

e d = (L-z?L z .)r)+ ]T L zf dy»Ar li + £ F^dy* A m + U^dy? A rj. 

Let £ e X P (Y) be an infinitesimal variational symmetry of the CR Cl,ti- Then, for 
all s E Sol(C) of the system if 

(1) 

d[(j\ S )yj d ^(z)(e)] = -[(L(j i s( X ))e+ J2 (L z ^(e))oj i s( X )+ 

(p,,i)eP 

+ (^ 1 ))°i 1 *(*'? = K« 1 )^) J 1 W)''i. (7-8) 

(2) If Jp(ir) — J 1 (7r) is the full 1-jet bundle, last equality takes the form 

-[{Ltfsix))? + {L z ^{e))o^s(x)], t = K(^)n M ) oj\s))*. 

(3) If, in addition, 11^ = 0, \i = 1, . . . ,m, then V( 6 j and for all solutions 
s G T(ir) the Noether conservation law holds: 

d[(j\ S )r J dL - n (z)(e)] = Kw'sixve + (i,^ 1 )) ^v*)]^ = o. 

Last statement of this theorem is the standard form of Noether Theorem in 
Lagrangian Field Theory. 

Theorem 10. Noether Theorem, RET, k = 0. Let C be a lifted covering 
constitutive relation of the RET type with 9^ = F^dy^ A r\i + U^dy^ A n, and let 
£ £ X(Y) be a variational symmetry of C. Then for all solutions s e r(7r) of the 
balance system ~k, 

7.3. Full invariance and the source charge. Let now condition (7.1) is fulfilled 
i.e. G is symmetry of both flux and modified source term = n M w M An+F^w^An 

of a covering constitutive relation C . Notice that this term depends only on the C 
but not on the lifting of the CR C to the covering constitutive relation. 

Then, for any £ e g C X(Jp(nj), lifting £ and the Poincare-Cartan form to 
Jp +1 (n) and using the fact that any (n+l)-form on J^ +1 {tt) is contact we have 

7rf +1) *<%e™ +1 = -i^+idQ^ 1 = -i^+iidn^Aoj^An+n^duj^An+diFpAuj 1 ; Arj+Ffa? At?)) = 

= -ip [(rf h n M + d v n„) A uj" A n + d h F^ A w? A n + d v Fl A wf Aij] = 

= -i^ [(diUf.dx 1 + dyll^) A ^ A n + d^dx 1 A A n + d v F^ A lo^ An] = 

= -i e A^Aj| + d v F^ A lo^ An] = -i e 2Con = Con. (7.9) 

Here, as before Con is a contact form while 2Con - two-contact form (see Sec. 2). 
We have used here the equalities dw M An — —dz^ A dx^ A rj = and the same for 
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efcjf A n. During this calculation we repeatedly used the equality dx^ A n = but 
we have not used the balance system. 

Reminding the decomposition d = dh + d v as the sum of vertical and horizontal 
differentials we see that the (n+l)-form 

Q(0 = i^Q n +1 = i e [\i^ An + FX A rj\ = i e K c (7.10) 

has contact differential and therefore, zero horizontal differential: dhQ(£) = 
0. Thus, it defines the class of cohomology [Q(£)] of the horizontal complex 

(AW0)' d *): 

[Q(0=^^c]ei^ o + r 1 (^ oo W)) 

If the class [Q] is zero, then Q(£) = + Con is the sum of horizontal differential 

of a form $(£) and a contact form and, therefore is also the sum of a differential 
d<f>(£) and the contact form 

Q(0 = + Con. 

Applying now the pullback by j k (s) (and j fc+1 (s) where appropriate) and recall- 
ing that we have not used the balance equation in the calculation above, we prove 
the following 

Theorem 11. Let, in addition to the conditions of Theorem 8, g is the Lie algebra 
of the variational (infinitesimal) symmetries of the source part of the constitutive 
relation, i.e. (7.2) is true. Then 

(1) 

dv6" +1 = Cont=> rfj fc (s)*z,-6" +1 = 
€ C w € C 

for all sections s of the configurational bundle n. 

(2) Form Q(£) = ir^ e+1 ^*i^kQ f ^' 1 = i^kKc defines the class of cohomology 
[Q(0] of the horizontal complex (/\* (J£° (it)) , d^) , 13 - g-charge of the 
source II = (IT^ 1 + F'wf) A n. 

(3) If the class [Q(Q] is zero for all (eg, there exists a form $c(£) linearly 
depending on £ (and its derivatives) such that Q(£) = d&c(£) + Con is the 
sum of a differential rf($c(0) an d the contact form. In this case locally 
(and in a topologically trivial domain, globally) 

/( S )*^ +1 = rfj fe+1 ( S )*$ c (0 

for this g* -valued n-form $c (that is natural to call the g-potential of the 
source Qq +1 ). 

Corollary 3. If G is the Lie group of symmetries of a constitutive relation C and if 
the class [Q] in the last Proposition is trivial, then there exists the (locally defined) 
g* -valued n-form <£>c : £ - *• A n (Jp (w)) such that in the conditions of Theorem 11 
the following conservation law will holds 

dijc(j k (s)(x)(o-j k + 1 (sy$ c (i)} = o 

for all solutions s £ r(7r) of the balance system Be and all £ G g. 
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7.4. Examples: Energy-Momentum Balance law, case of gauge symme- 
tries. 

Example 12. Energy-Momentum Balance Law, lifted CCR. Let v be a 

connection in the bundle 7r : Y — > X with the connection form K v = ® (dy^ — 

Then, the horizontal lift of the basic vector field d x i in Y is = d x % + T^d y ^. 
Its lift to the Lie vector field in Jp (jr) has the form 

Remark 14. Connection ^ is assumed to be compatible with the partial structure 
of Jp (tt) . In a case of if © if' -structure this means that the ^-horizontal lift £ to V 
of any vector field £ S X(X) preserving the i^©if'-structure should be such that its 
flow prolongation to Jp (w) is possible and preserves the (partial ) Cartan structure. 
In the case of an integrable AP-structure this requires that the components Tf of 
the connection v depends on the variables x 1 , dj 6 K and on j/ M (see [25], Sec. 8. 2) 
but not on the complemental variables in X. 

Consider a j/-homogeneous constitutive law C and the corresponding balance 
system Be- Let C be the lifted CCR corresponding to C and = F^uj^ Ar)i — 
H^dy^ A r\ — F^u)? A r\ - corresponding (corrected) Poincare-Cartan form. Calculate 
now 

% 6 >fi = i% F Z<^ ^Vk = F*(T? - - A 77 fcl . 

Next one has, 

d%6g = d k (F^-z^) V +F^-z^X G , kV -F^A Vki +Con = (d k +\ G , k )[F*(T%-z%)]rj+con, 

since j(s)*duj' 1 — dj*sLo^ 1 = 0. 

For the right side of balance law (7.7) we have 

Il^(d})r) + F^(e)v = [MIT - zft) + F*{d k T% - 

since = Tf - z£ . 

Applying the pullback by the j 2 s we get for the right side expression 

As a result, the balance law (7.7) takes, for the vector field the form 
(4+A G>fe )[F^r£-<)]o^ 

(7-11) 

Thus, the Energy-Momentum Tensor for the constitutive relation C has, thus, 
the form 

Ti=Fi(T?-z?). (7.12) 

In the left side of (7.11) stays the covariant divergence of the energy- momentum 
tensor. Right side of (7.11) represents the source/dissipation term. 
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Calculating derivative in the left side of (31.17) and canceling similar terms we 
reduce it to the equality 

[(d k + A G , fe )^](rf - zf )] o f s = U^jh)(Tf - z?)oji(s) 

& [(d k + \ G ,k)F* - nj(rf - zf )] o f s = o. (7.13) 

This gives the explicit presentation of the components of the energy-momentum 
balance law as the linear combination of the original balance laws with variable 
coefficients (rf — z%). 

Example 13. Energy-Momentum balance law, semi-P-Lagrangian case. 

Let Cl,yi be a semi-P-Lagrangian covering constitutive relation (see Theorem 3). 
Let us specialize the Noether balance law given in Theorem 8 above for the case of 
the vector fields £j = di + Pf <9 M introduced in the previous Example. This balance 
law takes, for a solution s : X — > Y of the balance system if, the form 

[(z4+ ]T Lj.(r?-z?)+ ]T n-W-z^&Wh = -(n^rf-^o/^ 

(7.14) 

In particular, for the case of the full 1-jet bundle J 1 (7r) and of zero connection 
rf = this balance law takes the conventional form 

[(L(j\s)(x))5i - L ^{j 1 (s)(x))s^ i )]- l = {U^){ 3 \s){x)). (7.15) 

Example 14. Pure gauge symmetry transformation. Let £ = ^d y v be a 
vertical (pure gauge) symmetry transformation of a constitutive relation C. Then, 
the Noether balance law corresponding to the vector field £ has the form 

[(eF;)ui(s))] ;xi = (epulis)) (7.16) 

- the secondary balance law defined by the vector field £ £ 5(C), see Lemma 9 
below. 

8. Secondary balance laws and the entropy principle. 

Let the system of PDE (-^r) be a balance system corresponding to a constitutive 
relation C. A natural question that generalizing the "entropy principle" of Con- 
tinuum Thermodynamics ( [2 1 j ) is - are there, except of the linear combinations 
of balance equations of the system *, nontrivial (see below) balance laws defined 
on the same bundle Jp{^) where the CR C is defined that are satisfied by all the 
solutions of the balance system 

Definition 10. Let C be a constitutive relation defined on a partial k-jet bundle 
Jp(7r) and let Be be corresponding balance system (if). We call a balance law 

(^°J?(s)), = 0°J? +1 (*)- (8-1) 

of the order q i.e. given by a (n+l) + (n+2)-form a = K l r)i + Qrj on the space 
Jpi 11 ) + ^p +1 ( 7r ) " ^ e secondary balance law for the system if if any solution 
s : X — > Y of the balance system (if) is at the same time solution of the balance 
law (8.1). 

Following is the list of four classes of secondary balance laws of the balance 
systems. 
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(1) An interesting class of the secondary balance laws, including the linear 
combinations of the balance laws of the system is determined by the 
following 

Lemma 9. Let a vertical vector field £ = £,^d y ^ £ V(w) is such that the 
condition FDiv{^) = F^d^ — is fulfilled. Then the balance law 

jI *( S )d(eFtvi) = eiv/ aeFt ) oj\ S )\ xi = (eu^oj\ S ) 

belongs to the space BCq- 

Proof. Follows from d(j$ *(s)^F?r)i) = £"d(# *(s)if »7i)+jp *(s)FDiv(£)r). 

□ 

(2) If a Lie group G is the (geometrical) symmetry group of the balance sys- 
tem if , it determines the family of the balance laws corresponding to the 
elements of Lie algebra g, see previous section. If the second order source 
part of these balance laws vanishes, we get the subspace BCc. s of the space 
BC C . 

(3) The entropy principle of Thermodynamics (see ([20] [21] [22]) requires 
that the entropy balance 

h*i = S, (8.2) 

with the entropy density h° , entropy flux h , A = 1,2,3, and entropy pro- 
duction plus source £ = E s + S p belongs to the space BCc of the balance 
system of a given theory. This requirement place a serious restrictions on 
the form of constitutive relation C and leads to the construction of a dual 
system in terms of Lagrange-Liu fields A M considered below (see [21] )■ Even 
more serious constitutional restriction is the proper II law of thermodynam- 
ics requiring that the entropy production E p is nonnegative, [2T| [22] . 

(4) If an integrable dynamical system can be formulated as a system of bal- 
ance equations (KdV equation is an example) then all the higher order 
conservation laws are the secondary balance laws of this balance system. 

8.1. Secondary balance laws for the RET balance systems. As an illus- 
tration of the geometrical approach to the balance systems and the example of 
application we present here some results from the paper [26) (in preparation). In 
the Rational Extended Thermodynamics developed by I.Mullcr, T.Ruggeri, and I. 
Shish-Liu, [21] the space of the dynamical fields is chosen to be large enough so 
that the constitutive relation C is defined on the bundle Y and, in our terms, is 
represented by the section of the bundle A^ +1 (Y) — > Y. 

dF° 

A constitutive relation C is called regular if the matrix is nondegenerate. 
For a regular CR C one can introduce new field variables w M = F® and to write the 
balance system in the form dtw 11 + d x iFf = IT^, ^ = 1, . . . , m. 

Theorem 12. 

Let C be a a constitutive relation of the RET type with F® = . 

(1) There is a bijection between 

(a) Functions h°{y^) satisfying to the (overdetermined) system of PDE 

(n^d y , - v?d y «)h° = 0, a + v, A = 1, 2, . . . ,n, (8.3) 
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where n A = F A y v(x,y)d yl * and such that the f vertical,) Hessian 



Hess(h°)(x ,y ) = 



d 2 h° 



lx=x ,y=y 



is nondegenerate in a neighborhood of a point (xo,yo), and 
(b) Secondary balance laws () with the functionally independent Lagrange- 
Liu multipliers A M = S-jj in a neighborhood of a point (xq, J/o)- 



(2) In terms of local vertical variables A M the density, 
corresponding secondary balance law are 



l,ux and the source of the 



'K°(x,y) = h (X(y)), 
K»(x,y) = ^(\(y)), v= 1,2,3, 



(8.4) 



8.2. Example: Cattaneo heat propagation law. Balance equations of this 
model have the form 



dt(pe) +div{q) = 0, 
d t (rq) + VA(tf) = -q. 



.5) 



Since p is not considered here as a dynamical variable, we merge it with the field 
e and from now on and till the end it will be omitted. On the other hand, in this 
model the the energy e depends on temperature $ and on the heat flux q (see [5] , 
Sec. 2. 1.2 or, by change of variables, temperature i9 = $(e,<?) will be considered as 
the function of dynamical variables. 

The secondary balance laws for Cattaneo model (including the original balance 
laws) have the form 



K 1 
K 2 



q 3 



5> 



5 3 



A° e - J* 



V 



v -i A J 

\°#e eq ds 



ffo\ 



+ 



rn 

V o y 

A°(i?y +K 1 (ti) 
\°($)q 2 + K 2 {$) 

Kj{\%\\q\\ 2 + k A ^)q A ) 



,0 



(#)q A }\ 



J 



(8.6) 



First and second balance laws in the system (8.6) are the balance laws of the original 
Cattaneo system. Third one one is the trivial balance law (see Sec. 2). 

Restrictions placed by the system (8.3) leads to the expression of internal energy 
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Forth column gives the balance law with the production term 



Aj(X°Aq\\ 2 +K A ^)q J 



A — 1 \0 



A 



2AS 



For a fixed d this expression may have constant sign as the function of q if and 
has the form 



only if all /^(i?) = 0. Therefore this is possible only if the internal energy (8.7) 



e = € eq (-d) + 



2A,* 2A ,, \A,*J 



(8.9) 



with some function A°($). As a result, Cattaneo model has the secondary balance 
law 



A u e 



\°^d S + -T(0)A^\°Jq\f 



AW 



-A^Ag||g[[ 3 . 



(8.10) 



with the production term that may have constant sign - nonnegative, 
provided 



A^Ag ^ 0. 



.11) 



Last inequality is the II law of thermodynamics for Cattaneo heat prop- 
agation model. 



9. Conclusion 

In this work we've presented a variational theory of system of balance equations, 
realization of the constitutive relations of such system as an abstract Legendre 
transformation, invariant form of the balance system and the Noether Theorem 
associating with the geometrical symmetries of the constitutive relation the corre- 
sponding balance law (and in the appropriate cases - the conservation law). Further 
development of this scheme including algebraic and geometrical structures related 
to a balance laws, some functorial properties of developed formalism, more detailed 
relations to the variational bicomplex will be presented elsewhere. Applications to 
the thermodynamical systems: 5-field thermodynamical system, 13-fields systems, 
description of secondary balance laws using the exterior differential systems (see 
[3]) will be presented in the forthcoming paper [25] . 



10. Appendix 1. Proof of Proposition 6. 

Proof. Let £ = £?di + + £^<9 z f be any vector field in J 1 (7r). Present any vector 
field in the space A^ 1 (J 1 (ir)) in the form 
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We have: 

C v Q n+1 = {di v + i^d)[(q^ + qldz?) A 77] = 
= d[(q„e + g* tf) A - fody" + g^f ) A CVi] + k* [(^ A dy» + dg* A dzf ) A 77] = 
= £"dg M A 77 + g^" A 77 + gdql A 77 + gj,d# A 77- 
- (dg p A dy» + dq^ A dz?) A ?m + fady* + A (fa + C^gM + 

+ i v [(dg M A dy» + dql A dz£) A 77] = 
= edq^ A r? + g^df" A 77 + £f dgj, A r? + g* d£<* A 77- 
- (dg M A dy" + dg^ A dz$) A ?m + ( 9i dy" + A di« G (D»?)]+ 

+ + A 77] + [(-£"dg M - tfdgj.) A 77] + [(dg„ A dy" + dg^ A dzf) A CVi] ■ 

(10.1) 

It is easy to sec that all terms containing dg M and dg^ will cancel and we get 

A* Q n+1 = A 77 + g^dtf A 77 + div G {l){q^ + q^dz?) A 77)]+ 

+ [(£**di/' 1 + ^dz< i ) A 77]. (10.2) 

Calculating differentials, leaving only vertical differential due to the presence of 
77 as common factor and equating result to zero we finally get 

C^Q n+1 = Iq^e + q^ + diVG&iq^ + qldz^ + i^dy^ + ^dz^Ar! = 

(10.3) 

Using the fact that the £ is projectable to X, so that £^ depend only on X and 
using relations 

^/ = f,y^ J + f^u) 
we reduce last relation to the following 

g4^+C^)+^(C^^ 

(10.4) 

From this it follows that the prolongation of £ to the A 2 l+1 (Jp (71-)) with the prop- 
erties listed in the Proposition exists, is unique and given by 

P" = -q l >&-ql% y ~-q v divG(Q, 

□ 
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